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, CHAPTER I 

INTRODUCTION 

I.l GENERAL 

The studies of electron momentiom distributions (EMD) in 
atoms, molecules and solids have received renewed attention 
during the last few years (see Epstein^ for a detailed recent 
review) . Many new techniques such as Y-scatterinq etc. have 
been developed extensively and employed to measure EMD* s expe- 
rimentally. A host of these measurements have since then been 
carried out on atoms, molecules and solids. A number of theore- 
tical calculations and interpretation of EMD' s have been repor- 

2 

ted. The first review on this subject was pxobli shed in 1933, 

1 3—7 

followed by a spate of review articles ' only after 1970. 

If the number of review articles is taken as a fair indication 



4 


of the activity in a given subject# one can say that studies of 
EMD* s in atoms, molecules and solids have been reactivated after 
overcoming a prolonged period of dormancy. 

1 3—7 

Since many exhaustive review articles ' on this subject 
are available, we shall deal in this chapter only with some 
salient features of experimental and theoretical studies of 
EMD' s and the most recent work done in this area. 

I . 2 ELECTRON MOMENTUM DISTRIBUTIONS AND COMPTON PROFILE 

The first experimental observation of the increase in wave- 
length of electromagnetic radiation on scattering by matter was 

s 

made by Eve with Y-rays. This was followed by similar studies 
9 10 11 

by Florance and Gray. Compton gave further experimental 

evidence of this phenomenon with X-rays, in the form of spectral 

plots and also a theoretical explanation of the scattering pro- 
1 2 

cess involved. An identical explanation was developed indepen- 

13 14 

dently by Debye. Compton showed through further experiments 

that the modified line was too broad to be explained by the non- 

15 

monochromatic nature of the incident- X-ray beam etc. Jauncey 

attributed this broadening to the fact that the electron is 

. 16 ''' 

initially not at rest but has a finite momentum. DuMond 
studied this Doppler broadening process in detail. 

When an electromagnetic radiation of wavelength x' is scatte- 
red inelastically by a stationary electron one obtains from the 
laws of conservation of energy and momentum a single modified 
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wavelength X' given by 


12,13 


X ' 


2h . 2 
— sin 
me 2 


( 1 . 1 ) 


where 0 is the scattering angle = 0.0243 A). However, when 

the dynamical nature of atomic or molecular electrons is taken 
into account, the scattered photon will carry some information 
about the EMD' s of the scatterer through Doppler broadening. 

In Fig. 1.1 the basic processes in Compton scattering are shown 
(vide ref. 3), taking into account dynamical nature of the elec- 
tron and assigning a linear momentxam p to it. Taking the momentum 
of the incident photon to be and that of the scattered photon 

to be hk' , from the law of conservation of momentum we have 


hk' - hk = p - p' 

'w-V W* 



( 1 . 2 ) 


and from the law of conservation of energy 


hclk'j-hc|k| ~ (p^-p* 2) 


i 


(1,3) 


Defining the scattering vector K = k' -k and assuming 
|id>- > I k| -|k*| , the final momentum p’ can be eliminated from 
(1,2) and one obtains • 


AX =X'-X 


“ sin^ (-^) - 
me 2 


2 X sin ( 2 ^ 


me 


(1.4) 




FfgJJ Schematic diagram ot the interaction of 

an X- fay with an electron of momentum p . 



Fig. 1.2 Doppler broadening of the Compton (he 
due to finite initiM momentum of 
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where is the projection of the electron’s initial momentuin 
along the scattering vector. This modified Compton wavelength 
and Doppler broadening of the Compton line due to finite ini- 
tial momentum of the electron are shown in Fig, 1,2. The inten- 
sity of the line corresponding to the component of momentum p 

z 

will be proportional to the probability of finding that parti- 
cular component of the momentum; 


J(p^) = / 0 (Px-Py-Pz* aPx'^Py 


(1.5) 


where P (p) is the momentum density. 

In the case of isotropic momentum distributions (1.5) 
reduces to 


CO 



where l(p) = / P (p) .p^ . sin 0 de d<) 

P P P 


( 1 . 6 ) 

(1.7) 


and q is the projection of the electron' s initial momentum pn 

the scattering vector. A quantum mechanical derivation of the 

17 

scattered intensity distribution was given by Kilby. Starting 

18 

from the Waller-Hartree theory of scattering; 


( 


da \ 

dw dn- 


(JL£) 


(~) 


S(K,«,) 


( 1 . 8 ) 
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where /jrj^ refers to the Thomson classical scattering cross 

section^ ^ and is given by; 


( 


da s 
dfi^Th 


- 2 2 


me 


(1.9) 


Here, e and e' refer to the polarization of incident and scatte- 
red photon respectively. and denote the angular frequencies 
of the initial and scattered photon, S(K, oj) is now expressed in 
terms of the quantum mechanical matrix elements, viz.. 


S(K,io ) = £1 

, 1 

1 


j I j 0>p6(E^-S^(-h w) '|\ 


r 

SA 


( 1 , 10) 


Here K is the scattering vector and i j> and i0> are the initial' 
and final states of the system. The ^-function ensures conserva~ 
tion of energy during the scattering process under consideration* 
If one makes the following assumptions: 


(i) scattering is by single electrons; 

(ii) the final state of the electron is a plane wave ; 

(iii) the energy transfer in the scattering process is 
considerably larger than Eg , the binding energy of 
electrons in the scatterer ; 

S(K,m) reduces to J(p„) given above (vide Eq, (1.5)) . 
These approximations form the basis of the so-called “impulse 
approximation". 
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1.3 IMPULSE APPROXIMATION ( IA) 

20 

The impulse approximation has been studied in detail. Eisen 
berger & Platzman^^^^^ calculated S(K,w) (in Eq. (1,8)) with 
exact continuum final state wavefunctions which involve 
confluent hypergeometric function for Is electron in hydrogenic 
atoms. They obtained an analytic expression for S(K,u) ^d 
showed that it reduces to that obtained with IA if 


exp { 



v3 t^y = 1 


(i.ii) 


By dropping of the commutator in (l.ll) one obtains an 
expression for S(K,<d) which is independent of V. According to 
Eisenberger and Platzman, the physics of the situation tells us 
that V has cancelled out from ihitial and final~state energies. 
Thus/ for short time of interaction the potential that the elec- 
tron is moving in is constant. They also showed that this 
approximation holds well for the atoms belonging to the first 
row of the periodic table. The correction to IA was shown by 

them to be ^( Eg /Ej^ ) ^ where Eg is the binding energy of the 

21 

electron and Eg / its recoil energy. Currat et al. and Bloch 

22 . 
and Mendelsohn have also studied the non-lA GP' s in more 

1 

detail. However, as pointed out by Epstein, IA retains its 
usefulness in the analysis of Y -ray CP data. The corrections 
to IA can be made, but they spoil the beautifully simple rela- 
tionship (Eqns (1.5) and (1.6)) to obtain CP from EMD' s. 
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In the next section we discuss briefly the experimental 
methods which have so far been employed to study EMD' s and CP' s 
for atoms / molecules and solids. 


1.4. EXPERIMENTAL DETERMINATION OF ELECTRON MOMEaSTTUM DISTRIBUTION 

There are six methods which have so far been employed in 
the experimental measurements of EMD' s of atoms# molecules and 
solids; 

a. X-ray Compton scattering. 

b. y -ray Compton scattering. 

c. High-Energy electron scattering. 

d. Positron annihilation angular correlation. 

e. (e, 2e) technique. 

f. Synchrotron radiation scattering. 

1 

Since a recent exhaustive review dealing with most of the 
above methods is available we shall discuss here briefly only 
the basic features of the above methods. Some of the recent 
experimental results are also included for the sake of complete- 
ness. 

1*4. A X-ray Compton Scattering 

A typical experimental set-up is shown in Fig. 1.3. Sources 
of radiation commonly used are MoKci (17.4 KeV) and AgKa (22.2 KeV) 
lines. A collimated X-ray beam is scattered by the sample and 
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then diffracted by an analyzing crystal to measure the wave- 
length (and hence energy) of the scattered radiation via Bragg's 
law of diffraction. The most severe restriction of X-ray Compton 
scattering is relatively low intensity of inelastically scattered 
radiation due to high photoelectric absorption cross-section. 

As a consequence, high statistical accuracy and good signal to 
noise ratios are not easily obtained. The incident radiation is 
non-monochromatic due to the doublet nature of X-rays from the 
sources discussed above and one cannot afford to lose the 
scattered intensity further by monochromatizing the incident 
beam. The CP has to be obtained point by point due to the use 
of crystal set for only one energy at a time. The background 
problem is quite complex in the case of X-rays and the resolu- 
tion is variable. X-rays cannot be effectively employed for 
elements beyond the first row due to high photoelectric absorp- 
tion. For higher elements the binding effects are also signifi- 
cant. A factor which favours X-rays over Y-rays is high inten- 

4 

sity of the incident radiation, typically 10 times that of 
Y“ray sources. X-rays were exclusively employed before the 
advent of Y-ray techniques for CP measurements (see Epstein^ and 
Cooper for reviews)-. In recent years reports on the use of X-ray 
technique to measure CP* s have virtually disappeared from the 
literature. ■ 

I.4.B Y -ray Scattering 

Development of solid-state detectors such as Ge(I.i) and 
Si(Li) has accelerated the growth of this technique. A typical 
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Y~ray CP set-up is shown in Fig. 1.4. In yray scattering 

experiment, one can measure all the points of CP simultaneously 

with a multichannel analyzer. The other chief advantage of 

Y-ray experiment is that the Compton to photoelectric cross- 

section ratio, which is proportional to (E/Z) is quite high 

and CP measurements for heavier elements can be carried out. 

Purity of the incident beam is also ensured automatically. The 

background problem is much easier to handle in the case of Y-rays 

than that for X-rays. A critical comparison of X-ray and -y-ray 

23 

Compton scattering is presented by Eisenberger and Reed. They 
have also given formulae for all corrections (involved in Y -ray 
CP measurements) including the relativisitc ones which become 

24 

significant for highly energetic Y“rays. Fukamachi and Hosoya 

have studied the optimum conditions in CP measurements by using 

a solid-state detector. They found that as the energy of the 

incident photon becomes higher, the optimum scattering angle 

o 

decreases from its value of 180 . However, the counting effici- 
ency of the detector decreases. They also estimated the effect 

25 

of both of these factors. More recently, Chaddah and Sahni 

have reported a careful study of several features governing the 

choice of Y -ray source, scattering angle, collimation etc. They 

have pointed out that optimum choice of scattering angle and 

collimation can offset the poorer penetration of low-energy 

241 

Y-rays. They also showed that between Am (59.54 KeV) and 

1 O '3 m 

■ Te (160 KeV) the former can be a worthwhile choice in 
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contrast to the earlier reports in the literature. 

A review of the earlier CP measurements employing Y-ray 

scattering is available. The Y“ray CP measurements carried out 

after publication of this review article include comparison of 

26 

CP's of isomer-pairs: methyl formate and acetic acid as well 

27 

as cyclopropane and propylene. An annular source of Y -rays 

28 

was employed in the measurement of CP of solid' KF. CP* s of a 

29 

host of ionic solids have been measured (see Paakhari et al. 

30 

for a review) . Sivaramu et al measured differential Compton 

scattering cross sections for different angles (from 3C? to 130°) 

by interposing cylinders of scatterersof different radii. Swamy 
313 2 

et al. * have obtained Compton spectrum due to L-shell elec- 
trons by employing a coincidence technique, Paakkari and 
3 3 

Merisalo have measured the CP of gaseous N 2 .’^ -ray CP's of LiP 

with scattering vector along <100>, < 110 > and < 111 > crystallo- 

34 35 

graphic axes have been reported. Manninen et al. have measur- 
ed CP' s of formamide and p-benzoquinone and compared them with 

36 

theoretical CP's. Reed et al. have obtained CP's of graphite 
and diamond with Y -ray scattering techniques and compared their 
results with theory. More recently CP of neopentane has been 
measured.” 

The Y-sources commonly used for CP measurements are 

59.54 keV ^^^Am and 160 keV With the latter source, 

23 

Eisenberger and Reed were able to measure CP of Krypton in 

38 

three days, A novel method of using a relatively short-lived 
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isotope of extremely high intensity has also been studied. 

198 

412 keV Au (half-life 2.7 days) of intensity 70Ciwas used 
for the measurements of CP’s of solid Al and Ge. 

1.4. C High-Energy Electron Scatterin g 

39—42 

Bonham and coworkers have developed this technique 

extensively. The energy- loss spectra of H 2 / He, Ne and Ar were 

39 

measured at low resolution with 45 keV electrons. Wellenstein 
40 

and Bonham replaced the solid-state detector with an electro- 
static differential velocity analyzer. The corrections to be 
applied for the CP measurements using this method are difficult 
to estimate and involve exchange as well as correlations between 
pairs of electronic momenta. However, this method has signifi- 
cant potential due to better resolution and high scattered inten- 

1 

si ties. Epstein has pointed out that it should be possible to 

obtain a statistical accuracy of better than 0,1% by this method 

for light elements. Details of experimental set-up are given lay 

41 42 

Wellenstein et al. ' CP’s of methane and ethylene have been 

43 

measured by Lahmam-Bennani et al . using high energy electron 

impact spectroscopy with 35 keV-electrons. Klapthor and Lee have 

44 , 

also measured CP of methane by this method. The CP’ s of ammonia 

45 

and carbon dioxide have also been recently measured. 

1 . 4 . D Positron Annihilation Angular Correlation 

De Benedetti et al. pointed out that EMD of positron- 
electron pairs can be estimated by measuring the angle between 
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two photons produced in a 2 T- annihilation of a positron and an 

electron. This method yields more detailed (three dimensional) 

information about the EMD‘ s than that obtainable with CP measure 

ments. Metal-ammonia solution and organic liquids were the 

systems of chemical interest studied earlier by this method. 

49 

Hogg et al. have reviewed the earlier work done with this 

50 

technique. A recent review which summarizes the positron anni 
hilation work done during 1972-74 is also available. 

1 . 4 . E (e, 2e) technicaie 

In this method, direction and kinetic energy of scattered 
as well as ejected electrons are measured in coincidence. By 
this method one obtains information about the binding energy as 
well as EMD of the scattering electron. Wei gold et al.^^ have 
attempted to develop this technique by using lower energy elec- 
trons and symmetric non-planar geometry. Epstein^ has pointed 
out that this method of obtaining EMD' s is yet to be developed 
but is attractive- in that the EMD' s of electrons in different 
orbitals may be measured separately by employing this technique. 

1 • 4 . F Synchrotron-Radi ation Scattering 

5 2 

Recently Cooper et al. have carried out exploratory 
Compton scattering studies on lithium fluoride, aluminium and 
copper employing monochromatize^ radiation from a synchrotron. 
They employed Ge(Li) detector to obtain energy spectrum of the 
scattered radiation. Typical time taken to accumulate '-'10, 000 
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counts at the Compton peak was '^'IS minutes in the case of 
lithium fluoride. However^ resolution obtained in this method 
is rather limited. At shorter wavelengths it is comparable to 
that obtained in the case of Y -ray CP measurements. Thus, this 
technique is suitable for carrying out CP measurements at high 
count rates and good signal to background ratio at modest resolu- 
tion. 

We feel that electron as well as synchrotron radiation 
scattering and (e, 2e) technique are poised for extensive develop- 
ment. However, y-rsy measurements would continue to be popular 
due to the simplicity of the set-up and a better understanding 
of the various corrections that are. to be applied to the experi- 
mental data. In the next section, we shall review various correc-* 
tions involved in Y -ray and X-ray Compton scattering. 

1.5 CORRECTIONS TO X-RAY AND Y-RAY COMPTON SCATTERING DATA 

A detailed account of various corrections to X-ray Compton 

3 23 

scattering is given by Cooper. Eisenberger and Reed have 
discussed various corrections applicable to Y-ray CP measurements. 
The corrections common to both methods are deconvolution and 
multiple scattering. The other corrections to be applied to 
X-ray CP data include those for: 

(i) wavelength dependence of absorption of X-rays in sample 
and analyzing crystal. 
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(ii) Wavelength dependence of polarization factor. 

(iii) Removal of white background. 

(iv) Doublet nature of the incident beam. 

(v) Divergence of the incident beam. 

(vi) Application of Breit-Dirac factor. 

Corrections applicable to Y-ray CP measurements are; 

(i) Background subtraction. 

(ii) Absorption of ‘^''-rays in the sample. 

(iii) Relativistic corrections. 

(iv) Energy dependence of the Compton cross section and 
other energy dependent factors. 

Deconvolution^ multiple scattering and relativistic correc 
tions have been studied extensively during the last few years. 

We shall give here a brief account of these three corrections, 

I *5. A Deconvo lution 

Any spectral measurement can be regarded mathematically 
as a convolution process 

h(x) =J g(x-z).f(z) dz (1.12) 

where h(x) is the measured spectrum, f(z) is the true spectrum 
and g(x-z) is the instrumental resolution function. Applying 
fourier transform (FT) techniques, one obtains 


F(C ) = H(c )/G(c ) 


(1.13) 



where P, H and G are FT' s of f, h and g respectively/ x and z 


refer to frequency domain and? to the time domain. Paatero 
53 

et al. have discussed three deconvolution schemes in terms of 
harmonic analysis. 


(a) Stokes" methods 

F(c) is computed as follows: 

F(?) = H(?)/G(c) for ? (1.14) 

F(? ) = O for ^ > 5 ^ 

f(x) is synthesized back from P(?) . 


(b) Generalized least-squares method: 


The problem here is to find f(t) such that 
/ ( h (x) - / g(x-t) f ( t) dt ) ^ dx + ^ f ( t)^ dt (1.15) 


has minimum value. The latter term ensures a smooth solution. 

(c) Successive approximations: 

The matrix equation corresponding to (1.12) is 
H =:GF (1.16) 


where H and F are the vectors containing measured result and 
the true spectrum, Q is a band matrix containing the instrumental 
resolution function. 

The exact solution of (1.16) is 
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which can be found iteratively by employing the scheme 


■'P'. 

'^1 


F 


i-1 


+ 


(H-GP._J 



(1.18) 


' 53 

Paatero et al. have given a detailed discussion of practical 
solution of the least-squares problem and stripping of tails in 
order that G in (1.15) does not become unwieldy to handle compu- 
tationally. 

The residual instrumentational function ( RIF) is the instru- 
mental function of the whole process, viz. measurement and decon- 
volution. Paatero et al. have discussed RIF for the above three 
methods of deconvolution. These RTF's show oscillatory character. 
However, it would be desirable to have an os dilation- free RIF. 
Another way to avoid these problems is to convolute the theore- 
tical data with the resolution function and compare directly with 

experiment after other corrections have been carried out. This 

26 

approach has been followed in recent measurements. 


I . 5 . B Multiple Scattering 

Of late multiple scattering (MS) corrections to the CP 

measurements have been discussed extensively. Neglecting the 

contribution due to MS introduces errors typically '-->10% in the 

peak height of the CP. The first attempt to evaluate the effect 

54 

of MS on CP was done by DuMond, who assumed the electron to foe 

55 

stationary in his treatment. Williams, Pattison and Cooper^^ ^ 

'■'56 

extended DuMond* s work to the case of moving electrons. Mdntire^ 
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calculated the MS contribution to beryllium CP for back-scatte- 
ring geometry. Pelsteiner et al.^'^ and Pattison et al.^® applied 

Monte Carlo technique to estimate effects of MS in water and 

-I . 59 

aluminium Compton scattering. Williams and Halonen have 

investigated qualitatively as well as quantitatively the ratio 

of double to single scattering. They also carried out Monte 

Carlo calculations to study this ratio quantitatively and found 

that the effective sample thickness is a useful parameter for 

6 C 

determining the contribution by double-scattering. Halonen etal. 
applied a Monte Carlo method to calculate spectral distribution 
of double scattering for aluminium with different sample thick- 
nesses. 

More recently Tanner and Epstein have studied effects 

OX 03 

in Compton scattering in a series of papers. They calcula- 

61 

ted the probabilities of single and multiple Compton scattering 

and angular distribution of scattered radiation analytically 

assuming Thomson scattering and specific sample geometry. They 
6 2 

also gave an analytic as well as numerical treatment of energy 
profiles for twice scattered photons from stationary electrons 
from a cylindrical sample. Total probabilities^ angular distri- 
bution and energy profiles of multiply scattered radiation were 

63 

evaluated by them by employing Monte Carlo techniques, 

59 

As pointed out by Williams et al. MS is the only major 
outstanding problem in CP measurements and considerable attention 
is being paid to study it extensively. 
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I • 5 . C Relativlsitc Corrections 

Relativistic corrections become significant in the Compton 
scattering of high-energy y-rays. Jau’'l:h and Rohrlich^"^ have 

■‘W 

given the expression for total Compton cross section within the 
relativistic framework. Eisenberger and Reed gave a relativis- 
tic treatment of the Compton cross section for the case of free 
electrons.' They also pointed out that a variety of relativis- 
tic formulae for the total Compton cross section have previously 

been employed and this would yield significantly different CP' s 

(b (S 

from a given experimental measurement. Ribberforts has presen- 
ted a relativistic treatment of differential cross section for 

Compton scattering against bound electron states for all scatter- 

S 1 

ring angles. He also studied this for anisotropic momentum 

S 8 

distributions and polarized photons. More recently he has 
obtained an expression for relativistic differential cross 
section for the case of a relativistic Maxwell-Boltzmann distri- 
bution of electrons. 

Having surveyed major experimental methods to measure EMD' s 
and some current developments in the studies of various correc- 
tions to the scattering data, we discuss theoretical calculations 
and interpretation of EMD' s and CP' s of atoms/ molecules and 
crystals in the following sections. 
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1.6 SOLUTION OF SCHRODINGER* s EQUATION IN MOMENTUM SPACE 
Starting from the Schrodinger' s equation 

^ Ti>(r) = (r) 


( 1 il9 ) 


where '3rt = + v(r) 


( 1 . 20 ) 

and multiplying both sides of (1.19) by exp(-i^.r) followed by 
integration over r, one obtains 


^X(p) +/e"^£*5v(r) T|j(r)dr= EX(p) 




( 1 . 21 ) 


where x(p) is the Fourier transform of 'P (r) . This can be 
formally written as 


^ x(p) V(p'-p) x(p' ) dp' = E x(p) 


( 1 . 22 ) 


where V(p'-p) = / V( r) . exp Pi (p' -p) , r~| dr 

i— *-«— L 


(1.23) 


(1.21) is an integral equation in momentum space. However for 
many-electron systems, the solution of (1.22) is a formidable 


69 


task and such calculations have so far been done only for helium*, 

+ 70 

H 2 and 

A simpler approach for obtaining EMD' s theoretically is to 
Fourier transform (FT) the co-ordinate space wave-functions. This 
approach has been used extensively (see Epstein^ for formulae for 
FT' s of Gaussian and Slater type orbitals & further references). 
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I • 7 CHEMICAL INTERPRETATION OF ELECTRON MOMENTUM DISTRIBUTIONS 
7 1 

Coulson did the pioneering work in the interpretation 
of molecular EMD' s and CP* s from the chemist' s point, of view. 

The EMD' s for diatomic molecules obtained from Heitler- London 
as well as MO wavefunctions were analyzed by him. He obtained 
the result that for both types of wavefunctions the electron is 
more likely to move perpendicular to the bond rather than along 
it. Coulson attributed this phenomenon to a "diffraction factor" 
which arises from the form of the wavefunction employed. Recently a 

5 

better physical picture has been given in terms of charge accu- 
mulation in the intemuclear region during bond formation. 

7 2 

Coulson also studied the EMD* s obtained from hybridized orbitals 
for hydrocarbons and found that increased p-orbital contribution 
in a sp^ hybrid orbital makes the momentum distribution broader. 
Epstein and Lipscomb and Epstein found that the directional 
characteristics of EMD* s are less pronounced in polyatomic mole- 
cules than in the diatomic cases due to significant overlap. 

74 

Epstein has examined the sigma-pi separability in momentum 

75 

space in the case of benzene. Henneker and Cade have attribu- 
ted the nearly spherical momentiam density in the case of LiF to 
the high ionic character of the bond in this diatomic molecule. 

Such a feature is absent in the case of covalent molecules. Tawil 
76 

and Langhoff analyzed the momentum densities in N 2 / 0 ^^ and H 2 CO 
and found that the density shifts from very low momentum values to 
intermediate momentum. Hbwever/ the densities for higher values 
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of the momentum retain essentially atomic nature. More recently 

77 

Kaijser and Lindner have analyzed the momentum densities for 

N 2 / CO and BF in terms of their sigma and pi symmetries for 

various states of neutral and singly ionized states. Kaijser 
7 8 

and Smith have calculated directional CP's for N 2 . They have 

examined the possibility of obtaining molecular bond lengths from 

directional Compton profile by studying the oscillatory nature of 

79 

the Compton profile along the bond. Ahlberg and Lindner have 
examined the correlation between electrons with parallel spins 
(Fermi correlation) by introducing a correlation coefficient. 

They found that the same rules govern the Fermi correlation 
between electrons in position and momentum space. 

80 

Pattison et al. have described a new method of interpre- 
ting CP data using the autocorrelation function B(r) of the one- 

Vi»*\ 

electron wavef unction 'J' (r) 

B(r) = /'P(R) .'P (r+ R) dR (1.24) 

and showed that B(r) can be obtained in any particular direction 

by taking one-dimensional FT of the corresponding directional CP. 

They applied this model to LiP and showed that ionic radii can be 

81 

determined by this approach. More recently Schulke has studied 
the possibility of mapping the Fermi surface from one- dimensional 

FT of CP' s. He has also given theorems on one dimensional FT of 

' 82 ■ ' ' ' ' 

CP's and their applications. 



26 


1 . 8 MOMENTUM EXPECTATION VALUES 

83 '4 8^ 

Coulson, Benesch and Smith and Epstein have indepen- 
dently derived the relationship between various < p^ > expectation 
Values and CP. <p^> is given by 

00 

< > = 2in+l)J q^J(q) dq (1.25) 

84 

Epstein speculated that it may be possible to correlate < p > 

with some molecular property such as shielding in nuclear magne- 

2 

tic resonance. < p > is the kinetic energy and relativistic co- 

4 

rrections are proportional to <p > . Epstein also studied the 

possibility of calculating these momentum expectation values 

8 3 

from CP data. Coulson showed that the energy of a molecular 
system can be written by employing virial theorem, as 

®mol " 'V 

He pointed out that if individual-component CP' s are additive 
and their sum yields the total molecular CP(as in Localised 
Molecular Orbital (LMC) approach) is also additive. Thus 

the additivity of bond energies used extensively in classical 
thermochemistry is justified in the LMO approach. 

1.9 LOCALIS ED-MOLECULAR ORBITAL STUDIES OF COMPTON PROFILES 

85 ■ ■ 

Hicks was the first to visualize the use of LMO* s in 

CP calculations. He stated that it may be possible to predict 
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CP's of a large number of molecles from CP's for relatively 

7 4 

smaller number of molecules. Epstein obtained li-lO contribu- 
tions for a large number of bonds and lone pairs from Edmiston- 

86 

Ruedenberg- localized MO's. However he used an average LMO 

approach while dealing with multiple bonds. He also showed that 

the CP' s obtained by adding up the individual LMO contributions 

agree well with those obtained from a full calculation. 

87 

Eisenberger and Marra utilized this concept of LMO' s and 
derived the LMO profiles for C-H, C-C and C=C bonds from experi- 
mental X-ray CP data for a large ntamber of hydrocarbons. More 

88 89 

recently Smith and Whangbo ' have also carried out theoreti- 
cal calculations of LMO contributions to CP' s but assuming diffe- 
rent contributions for single and multiply bonded cases such as 
C-C and C=C. This approach can be carried one step further by 
considering the state of hybridization of each atom in a bond. 


1. 10 OTHER RECENT C/^J^CCLATIONS OF COMPTON PROFILES OF ATOMS , 
MOLECULES AND CRYSTALS 

A thorough review of theoretical calculations of CP' s is 
given by Epstein.^* This review covers generally the work done 
upto 1973. Therefore we shall give here an account of only the 
more recent work done in this area. 

90 

Sabin and Trickey ■ have calculated CP's from HFS atomic 

wavefunctions for various values of the Slater exchange parameter, 

91 

a . Biggs, Mendelsohn and Mann have prepared extensive tables 
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of Hartree-Fock CP's for elements upto atomic number 102 for 

numerical wavefunctions . They have calculated non-relativistic 

CP' s for atomic numbers Z < 36 and relativistic CP's for higher 

9 2 

elements. More recently Benesch has computed CP's for neutral 
atoms As (Z= 33) through Yb (Z = 70) . These tables^^"^^ will 
serve as reference work for atomic CP* s. 

93 94 95-97 

Smith and Brown ' and Benesch and Smith have inves- 

tigated the role of correlation in EMD' s and CP' s using a density 

-r 94 

matrix approach for Li, Li , Be, H 2 and Ne, Smith and Brown 

found that correlation plays a relatively minor role in EMD' s of 

closed-shell systems such as iielium and neon. Correlation, in 

general/ increases J(o) and has little effect on the tail of the 
98 

CP, Benesch has computed EMD' s and CP' s from two-electron 
atomic wavefunctions containing exponential correlation terras. 

He applied Fourier integral techniques to generate non-diagonal 
first order density matrices. He also pointed out that the 
techniques used in his analysis are general and can be applied 
to larger systems. However, lengthy algebra limits their use- 
fulness. 

A review of experimental and theoretical work done in CP's 

29 

of ionic crystals has been given by Paakkari et al. More 

99 100 

recently, Aikala / has presented a detailed account of theore- 
tical calculations of directional and spherically averaged CP' s 

using a linear combination of atomic orbitals (LCAO) based oh 

101 102 

Lowdin' s orthonormalization procedure. Seth et al. have 
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calculated LCAO band structure and CP of TiC. Seth and Ellis^*^^ 
have computed EMD' s and CP's for diamond, silicon and- silicon 
carbide with approximate HFS crystal wavefunctions. 

104 

Thulstrump et al . have given formulae for calculation 

of anisotropic CP's from LCAO-type molecular orbitals. Ramirez 
105 

et al. have theoretically calculated CP anisotropies in mole- 

106 

cules and solids, Snyder and Weber have computed CP of water 

from a SCF-MO wavefunction in a double-zeta Gaussian basis-set. . - 
107 

Whangbo et al. have estimated the effect of hydrogen bpnding 

108 

on CP of water. Smith et al. have investigated the roie of • 

correlation on the CP of water. They found that the: contribution 

of electron correlation to CP is not insignificant, but minor in 

comparison to the effects of basis-set extension. Ahlenius and 

Lindner^ have calculated CP's and<p^?- values for water, ammonia 

and methane using large-scale DZ wavefxanctions . Ahlenius and 

Lindner^ have calculated valence CP's from semi-empirical wave- 

functions such as CNDO/2 , INDO, EHT etc. and discussed the bond 

111-113 

profile transferability. Hirst and Liebmann have calcu- 

lated CP's for 18-, 24- and 42-electron series using MO wave- 

114 

functions employing Gaussian basis-sets, Ulsh et al. as well 

115 

as Braun-Keller and Epstein have analyzed vibrational correc- 

1 16 

tions to the EMD* s in the case of using Wang’s wavefunction. 
Ulsh et al. found that applying the vibrational correction 
increases the discrepancy betwe^ the experiment and theory. 
Braim-Keller and Epstein, on the other hand, pointed out that 
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the sign of the vibrational correction is very sensitive to 
small amount of variation of the electronic wavef'unctioh used 
in their study. 

I. 11 SCOPE OF THE PRESENT WORK 

It can be seen from the review of the studies of EMD* s in 

atoms, molecules and solids presented in Sections I.l through 

I . 10 that a vast amount of work has been carried out in this 

area. However, very few studies have been carried out on atoms 

117 

using the Slater exchange potential. It is also possible to 

simulate the wave functions of ions in crystals using a spheri- 
cal potential around the ion. In Chapter II, CP's of crystal- 
ions obtained from Hartree-Fock and Hartree-Fock-Slater wave- 

functions will be discussed. A spherical potential first proposed 
118 

by Watson has been employed in this model to simulate the 
behaviour of ions in crystals. 

The present work attempts to examine EMD' s and CP' s obtained 
from models typically employed in chemistry for studying molecules. 
EMD' s from two-centre VB wavef-unctions including ionic terms will 
be analyzed in Chapter HI. It may be noted that the only other 

71 

theoretical analysis of VB wavefunction was carried out in 1941. 

In Chapter IV, EMD' s and CP's from floating spherical Gaussian 

119 

orbital (PSGO) model of molecular structure, which has been 

extremely popular ab-initio model among chemists, will be^^^^ ^ ^ 
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discussed. It may be pointed out here that FSGO’ s have not so 
far been used for the calculation of EMD' s. In Chapter V, CP's 
of a few hydrocarbons computed from semi- empirical _Complete 
Neglect of _pif f erential Overlap ( CNDO) wavefunctions will 
be examined with a view to study bond- additivity of CP's. As 
noted in Section 1.9/ UMO approach to CP's has been recently 

shown to be promising. In Chapter VI, CP' s of acetone and allyl 

■ ’ 241 

alcohol obtained by the author by scattering of 59.54 keV Am 

Y-rays will be presented end results compared with LMO-theoreti- 

7 4 

cal -CP's. A method to evaluate atomic and molecular energies 
from experimental CP's by fitting the CP data to: a linear combi- 
nation of Gaussians will.be presented in Chapter VII. The 
energies thus obtained from experimental CP' s will then be 
compared with the experimental and HP theoretical ones. The 
possibility of employing this method to obtain accurate molecu- 
lar energies will then be discussed. We have recently observed 
a remarkable empirical correlation between the CP and molecular 
energy in iso-electronic series. Results of this empirical 
correlation and their implications will be discussed in 
Chapter VIII. 

1.12 SUMMARY AND CONCLUDING -REMARKS 

In the present chapter, we have given a brief account of 
phenomenon of Compton scattering for non- station ary electrons. 

An account of impulse approximation and various assiimptions 
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underlying it was presented in Section 1.3. Various experi- 
mental methods of determining s of atoms/ molecules and 

solids were surveyed in Section 1.4, followed by recent work 
done on various corrections involved in the CP measurements. 
Direct solution of Schrodinger' s equation in momentum space 
was briefly discussed in 1.6. An . account of chemical interpre- 
tation of EMD' s and momentum expectation-values was given in 
the next two sections. In 1.9, the IMO approach to molecular 
CP's was discussed, followed by a review in I. 10 of recent compu- 
tations of SMD' s of atoms, molecules and cirystals. Scope of the 
work presented in this thesis was discussed in Section 1. 11,. 

Recent theoretical and experimental work in the area of 
EMD* s and CP' s clearly point to the importance of these re- 
emerging topics . Further work may bring us closer to the 
answer for the classic question posed by Mulliken several years, 
ago: "What are the electrons doing in molecules?" 
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CHAPTER II 

COMPTON PROFILES OF FREE AND CRYSTAL- IONS 
WITH HARTREE-FOCK AND HARTREE- 
FOCK-SLATER WAVEPUNCTIONS 

II. 1 INTRODUCTION 

Recently in atomic structure theory^ models in which the 
potential energy jl,s a local functional only of the density,p ^ 
have been studied extensively (see Slater^ for an extensive 
review). The p ^ approximation introduced by Slater is pro- 
bably the most popular one. Writing the Hamiltonian operator in 
the case of atomic systems as 

T + (2.1) 

where T, and represent the kinetic energy. Coulomb energy 

and exchange energy operators respectively. The exchange opera- 
tor within the Slater approximation expressed in Rydberg units 
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i s given by 


V. 


X 


3a 

2 


IT 


(r)} 


1/3 


( 2 . 2 ) 


Here a is the Slater exchange parameter. The value of o=l has 
been proposed by Slater while the value ofa =2/3 has been pro- 
posed by Caspar^ as well as Kohn and Sham"^. From model calcula- 
tions on atomic systems and by comparison with Hartree-Fock cal- 
culations the value a=i has been shown to overestimate the 
exchange potential whereas the value of a =2/3 mderestimates it. 
Several other recipes for choosing the value of a have been given 

Lindgren and Schwarz^ suggest the use of different a's for 

6 

different shells, whereas Schwarz recommend that value of “ for 
which the energy from the Hartree-Fock Slater (HFS) wavef unction 
agrees closely with the corresponding Hartree-Fock (HF) energy. 
Several properties of atoms , molecules and crystals have been 
studied with these models which employ this local fionctional 
of density. However, very few attempts have so far been report- 
ed in the literature regarding the calculation of electron momen- 
tum distributions (EMD) and Compton profiles (CP) of atomic 
systems employing Slater exchange potential. Singh and 

Smith"^ have studied the electronic radial momentum distributions 
of two-, three- and four-electron ions using HF and HFS wavefunc- 

tions. They found that the HFS momentum distributions depend 

3 ' ' ' ' ' ' ' 

significantly on the value of a used. Euwema and Surrat 
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compared the HP and HFS CP' s for Ne, Ar and Kr. They also 
found that there is a strong dependence of CP's on the value 
of dused. These authors also found that for Ar and Kr when 
the value of. a is chosen such that the virial theorem is satis- 
fied (a =a ), fairly good agreement between the corresponding 

HFS and HP CP' s results. However# for Ne, this is not the case- 

9 

More recently, Sabin and Trickey have compared the atomic CP' s 
of He, Li, Be, B and Ne obtained from four model local den- 
sity functionals. They found that the CP employing the model 
for which “"“yip gives smallest standard deviation from the CP 
obtained from correlated wavefunctions for the five atoms studied. 
However, to the best of our knowledge no extensive calculations 
of CP's of free ions using HFS wavefunctions have been reported 
in the literature so far. 

On the experimental side, CP' s of several ionic solids 
have recently been measured (see Aikala^^ for a review) .An 
overlap model^*^ employing crystals constructed by means of 
Lowdin' s symmetrical orthonormalization method^^ has been used 
to calculate CP's of these ionic crystals. This model improves upon 
the free ion results and gives better agreement between the theo- 
retical and experimental CP's. Ihere is yet another simplified 

model, which simulates the behaviour of ions in crystals, due to 
12 

Watson. In this model, a spherically symmetric external poten- 
tial around the ion is employed to simulate the environment of 
the ion in a crystal. This additional potential is of the form 
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V = n/r^ 

for 

^1^0 

{2.3a) 

= n/r 

for 

r > r 

<2. 3b) 


o 


is the radius of the sphere surrounding the ion and the 

number n of charges is chosen as n , = n. , where n. is 

the charge on the ion. This potential is shovro schematically 

in Pig, 2,1. Watson model employing the spherical potential 

13 

has been used by Paschalis and Weiss to calculate various 

properties of crystals such as diamagnetic susceptibility, 

dipole polarizability etc. These authors employed STO basis set 

within the Hartree-Fock framework in their calculation. More 

14 

recently Sen and Narasimhan have studied the Stemheimer anti- 
shielding factors and dipole polarizabilities of ions using 
crystal-ion HFS wavefunctions . Thus, it was felt that it would 
be worthwhile studying the CP' s obtained with the Watson model 
using HF and HFS wavefunctions and compare the CP's of ionic 
crystals thus obtained with other theoretical and/or experimen- 
tal data. When the present work was in progress, a report on 

15 

measurement of CP of solid KP was published by Weyrich. He 

also employed the crystal ion wavefunctions of Paschalis end 

13 ' ' ' " ■ ^ 

Weiss for calculating crystal ion CP' s of K and F and found 

that the crystal CP thus obtained agreed better with the experi- 
mental one than the corresponding free ion CP. 

We had independently confirmed this observation of Weyrich 
in the case of CP' s of several other ionic crystals calculated 
from free and crystal ion wavefmctions of Paschalis and Weiss t 
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+ 2 ^ 

In the following sections/ we present the CP' s of Li , Be , 

2+ + -f — 

Mg , Na , K and Cl calculated from HF and HFS wavefunctions 

— 2 - 

for both the free and crystal-ion cases. CP' s of F and 0 
using HF free as well as crystal-ion wavef motions and CP of free 
Br within the HFS framework are also presented. The CP' s of 
ionic crystals thus obtained are then compared with the corres- 
ponding experimental CP' s. The variation of HFS ionic CP* s with 
the Slater exchange parameter a is also studied. 


II. 2 METHODS OF CALCULATION 

The HF free as well as crystal ion wavefunctions were 

13 

taken from the tables of Paschalis and Weiss. The sphere 

16 

radius r^ closest to the Pauling ionic radius was chosen. Radial 
momentum densities I(p) were computed analytically using the 

formulae given by Epstein . CP' s, J ( q) , within the impulse approxi- 

1 8 * 

mation/ were calculated numerically from these momentum densi- 
ties by nxamerical integration using Simpson' s rules 

00 

J{q) = ^ J l(p) dp /p (2.4) 

N 

The HFS free- ion wavefunctions were generated in the 

19 

present calculations with a modified version of Herman-Skillman 
atomic structure program. The crystal-ion wavefunctions were 
generated by using the crystal-potential in equations ( 2.3 a &b) , 

20 

The HFS programs developed in 0*^^^ laboratory by Drs. K.M. S. Saxma 
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21 

and K.D, Sen were usecS in these calculations. A 441-point 
mesh was employed to carry out the integration in the equation 

I (p) = E !/ r^ r) j (pr) drj ^.4trP^ (2.5) 

nl ^1 

Here j (pr) is the Bessel function and a , is the HFS orbital - 
The general subroutine developed by Dr. Saxena was used to carry 
out this integration numerically. The CP's were calculated 
from I(p) ( which was obtained with equation (2,5))by employing 
equation (2.4). The CP' s were normalized to the corresponding 
number of electrons in the range O^q^lS. Calculation of CP 
of a shell took a computer time of 5-8- minutes on our IBM 7044/ 
1401 system, 

II. 3 RESULTS AND DISCUSSION 

19 

It is well-known that for® =1 and with Latter-Correction 
to give proper asymptoti c behaviour, the HFS wavefunctions agree 
fairly well with the corresponding HF wavefunctions. When Latter 
correction is employed along with the value of o=2/3 it is 
known that poorer HFS wavefunctions result. We may also note 
here that the free— ion ILFS wavefunctions and their a- dependence 
have already been studied extensively in the literature.^ However, 
no study of the variation of a in the HFS framework with the 
Watson sphere potential for ions has been reported so far in 
the literature. Hence, in the present program of work we have 
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included a study of the variation of a and its effect on crystal 

4 * 

ion HFS wavef unctions. Weyrich has noted in the case of K and 

F ions that the influence of the radius r of the Watson sphere 

o 

is small compared to the main effect of the charged Watson 

sphere. Hence, in all the ions studied here, we have fixed the 

16 

r^ value, at the corresponding Pauling ionic-radius and not 
examined the effect of variation of r^ on the crystal-ion CP's. 

Tables 2.1 to 2.5 present our results on cations while 
Tables 2.6 to 2,9 present our results on anions. Table 2.1 
shows the CP's of Li calculated from HF as well as HFS (a = l) 
free-ion wavefunctions. The HF and HFS (a=l, 0.77, and 2/3) 
crystal-ion CP's are also shown in this Table. Tables 2.2, 2.3, 

2.4, 2.5 and 2.8 show the CP's of Be^"^, Na"^, Mg^”^, and Cl" 
ions respectively calculated from HF and HFS (a=l) free ion 
wavefunctions. These tables also show the CP's of the above 
crystal-ions with HF as well as HFS (with a -variation) wavefunc- 
tions. It is known that the convergence in HFS calculations is 

2 - 

slow in the case of F and 0 ions. Hence these ions were not 
studied within the HFS framework. In Table 2,9 CP of Br calcu- 
lated from HFS free-ion wavefunction is presented. 

We note that the agreement between the free-ion CP's cal- 
culated with HF and HFS (a=l), wavef met ions is fairly good in 

4 - 2 ^ 4 * 24 * 4 “ — 

all the cases studied, viz., Li , Be , Na , Mg , K and Cl . 

The CP' s of crystals LiF, NaF, KF, Lid, KCl, NaCl, LiBr, NaBr, 

KBr, BeG and MgO calculated by using these free-ion CP' s are 
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Table 2.1 




Compton 

Profile 

of Li”^ 




HF 

HFS 

HP 

HFS 

HFS 

HFS 

Q 

free 

free 

r^=1.13 

r =1.13 
o 

a= 1 

r =1.13 
o 

a =0.77 

r =1.13 
o 

a =2/3 

0.0 

0.656 

0.664 

0.674 

0.716 

0.738 

0.750 

0.1 

0.652 

0.660 

0.672 

0.712 

0.734 

0.746 

0.2 

0.644 

0.650 

0.660 

0.700 

0.722 

0.732 

0,3 

0.628 

0.636 

0.644 

0.680 

0.700 

0.710 

0.4 

0.608 

0-614 

0.6 20 

0.656 

0.672 

0.680 

0.5 . 

0.584 

0.580 

0.594 

0.624 

0.640 

0.644 

0.6 

0.556 

0.560 

0.562 

0.590 

0.600 

0.606 

0.7 

0.526 

0.5 28 

0.528 

0.552 

0.560 

0.56 2 

0.8 

0.494 

0.496 

0.494 

0.512 

0.518 

0.5 20 

0.9 

0.460 

0.460 

0.458 

0.472 

0.476 

0.476 

1,0 

0.426 

0.426 

0.422 

0.432 

0.432 

0.432 

1.2 

0.360 

0.359 

0.354 

0.354 

0.352 

0.358 

1.4 

0.298 

0.296 

0.292 

0 . 286 

0.282 

0.280 

1 .6 

0.244 

0.242 

0.23 8 

0.228 

0.222 

0.2 20 

1.8 

0.198 

0.196 

0.194 

0.180 

0.174 

0.17 2 

2.0 

0.160 

0.156 

0.156 

0.142 

0.136 

0.134 

3 .0 

0.054 

0.05 2 

0.05 2 

0.044 

0.040 

0.038 

4 .0 

0.018 

0.020 

0.018 

0.014 

0.014 

0.013 

5 .O 

0.008 

0.008 

0.008 

0.006 

0.006 

0.005 
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Table 2,2 


2 + 

Compton Profile of Be 


q 

HP 

free 

HPS 

free 

HP 

r =1.11 
o 

HPS 

r =1.11 
o 

a =1 

HPS 

r =1.11 
o 

a =0.77 

0.0 

0.472 

0.476 

0.5 20 

0.583 

0.590 

0.1 

0.472 

0.476 

0.518 

0.581 

0.587 

0.2 

0.468 

0.472 

0.514 

0.577 

0.581 

0.3 

0.462 

0.466 

0.504 

0.565 

0.571 

0.4 

0.456 

0.458 

0.492 

0.551 

0.556 

0.5 

0.446 

0,449 

0.478 

0.534 

0.538 

0.6 

0.434 

0.437 

0.462 

0.513 

0.518 

0.7 

0.422 

0.424 

0.444 

0.491 

0,495 

0,8 

0.408 

0.410 

0.424 

0.467 

0.470 

0.9 

0.392 

0.393 

0.402 

0,442 

0.444 

1.0 

0.376 

0.376 

0.380 

0.415 

0.417 

1.2 

0,342 

0.343 

0.338 

0.362 

0.362 

1.4 

0.306 

0,306 

0.296 

0.310 

0.312 

1.6 

0 . 27 2 

0.271 

0.258 

0.26 2 

0.262 

1 .8 

0.238 

0.237 

0.222 

0.218 

0.217 

2.0 

0 . 208 

0.206 

0.192 

0.181 

0.180 

3.0 

0.096 

0.095 

0.088 

0.066 

0,065 

4.0 

0.044 

0.042 

0.040 

0.025 

0.024 

5 .0 

0.020 

0.020 

0.020 

0.010 

0.010 
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'I'cOdIg 2#3 

+ 

Compton Profile of Na 


HF HFS HF HFS HFS 

free free r =1.80 r =1.80 r =1.80 

o o o 

a =1 a =0.73 


0.0 2.285 2.236 


0.1 

2.282 

2.236 

0.2 

2.268 

2,222 

0.3 

2.240 

2.197 

0.4 

2.210 

2.163 

0.5 

2.118 

2.118 

0.6 

2.104 

2.063 

0.7 

2.040 

1.999 

0.8 

1.962 

1.925 

0.9 

1.882 

1.842 

1.0 

1.792 

1.755 

1.2 

1.602 

1.572 

1.4 

1.404 

1.385 

1.6 

1.218 

1.208 

1.8 

1.048 

1.048 

2.0 

0.896 

0.901 

3 .0 

0.424 

0.443 

4.0 

0.23 2 

0.246 

5.0 

0.146 

0.152 


2.294 

2.328 

2.509 

2.288 

2.3 24 

2.503 

2.274 

2.310 

2.485 

2.246 

2.288 

2.448 

2.216 

2.250 

2.395 

2.124 

2.201 

2.330 

2.110 

2.135 

2.246 

2.040 

2.065 

2.149 

1.961 

1.981 

2.041 

1.882 

1.884 

1.927 

1.792 

1.787 

1.809 

1.602 

1.575 

1.570 

1.404 

1.367 

1.341 

1.212 

1.178 

1.132 

1.042 

1.018 

0.952 

0.890 

0.876 

0 . 804 

0.422 

0.419 

0.388 

0 . 23 2 

0.240 

0.224 

0.146 

0.150 

0.142 


li.T. Kmnm 

^ A Mil’S 
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Table 2.4 

2 

Compton Profile of Mg 


q 

HF 

free 

HFS 

free 

a=l 

HF 

r =1 . 23 
o 

HFS 

r =1.23 
o 

HFS 

r =1.23 
o 

« =0.73 

HFS 

r =1.23 
o 

a= 2/3- 

0.0 

1.974 

1.941 

2.002 

2.182 

2.289 

2.315 

0.1 

1.970 

1.937 

2.000 

2.178 

2.285 

2.311 

0.2 

1.960 

1.927 

1.990 

2.165 

' 2.270 

2.295 

0.3 

1.946 

1.912 

1.976 

2.142 

2.245 

2. 267 

0,4 

1.924 

1.890 

1.948 

2.108 

2.209 

2.230 

0.5 

1.884 

1.860 

1.918 

2.067 

2.161 

2.182 

0.6 

1.858 

1.822 

1.882 

2.016 

2.102 

2.121 

0.7 

1.816 

1.784 

1.840 

1.955 

2.033 

2.050 

0.8 

1.770 

1.736 

1.788 

1.866 

1.953 

1.968 

0.9 

1.714 

1.680 

1.7 26 

1.809 

1.865 

1.877 

1 .0 

1.650 

1.624 

1.668 

1.7 25 

1.770 

1.779 

1.2 

1.522 

1.497 

1.534 

1.547 

1.569 

1.57 2 

1 .4 

1.384 

1.364 

1.384 

1.366 

1.367 

1.365 

1 .6 

1.242 

1.226 

1.236 

1.191 

1,176 

1.171 

1 .8 

1 . 104 

1.093 

1.092 

1.030 

1.005 

0.997 

2.0 

0.97 2 

0.968 

0.960 

0.887 

0.858 

0.847 

3,0 

0.278 

0.290 

0.266 

0.253 

0.222 

0.233 

4.0 

0.170 

0,182 

0.170 

0.164 

0.159 

0.155 

5.0 

0.116 

0.122 

0.116 

0,110 

0.107 

0.104 
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Table 2.5 





Compton 

Profile 

i* 

of K 



q 

HF 

free 

HFS 

free 

a=i 

HF 

r =2.51 
o 

HFS 

r =2.51 
o 

a= 2/3 

HFS 

r =2.51 
o 

a =0.71 

HFS 

r =2.51 
o 

a =1 

HFS 

r =2.51 
o 

a =1.1 

0.0 

4.428 

4.477 

4.456 

4.857 

4.805 

4.534 

4.468 

0.1 

4.406 

4.459 

4.438 

4.856 

4.786 

4.508 

4.452 

0.2 

4.358 

4.416 

4.387 

4.769 

4.719 

4.470 

4.400 

0.3 

4.282 

4.339 

4.295 

4.646 

4.601 

4.374 

4.308 

0.4 

4.143 

4.214 

4.168 

4.462 

4.424 

4.225 

4.182 

0.5 

3 .978 

4.035 

3.996 

■ 4.216 

4.188 

4.043 

3.993 

0.6 

3.774 

3 .802 

3.789 

3.921 

3.903 

3.796 

3.774 

0.7 

3.541 

3.542 

3.546 

3.594 

3.588 

3.537 

3.527 

0 .8 

3.288 

3.273 

3.282 

2.254 

3.257 

3.376 

3.253 

0.9 

3.024 

2.984 

3.018 

2.926 

2.937 

2.970 

2.991 

1 .0 

2.762 

2.706 

2.748 

2.617 

2.633 

2.701 

2.7 27 

1.2 

2.275 

2.213 

2.251 

2.093 

2.114 

2.213 

2.146 

1.4 

1.864 

1.820 

1.852 

1.701 

1.721 

1.820 

1.854 

1 ,6 

1.550 

1.5 23 

1.540 

1.424 

1.440 

1.5 23 

1.555 

1.8 

1.316 

1.309 

1.311 

1.242 

1.253 

1.309 

1.331 

2.0 

1.146 

1.141 

1.145 

1.081 

1.087 

1.141 

1.16 2 

3.0 

0.747 

0.756 

0.746 

0.740 

0.745 

0.756 

0.757 

4.0 

0.541 

0.542 

0.540 

0.537 

0.540 

0.547 

0.548 

5.0 

0.391 

0.386 

0.390 

0.384 

0.385 

0.386 

0.389 
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Table 2.6 




Compton 

Profile 

of F 



q 

HF 

free 

HP 

r =2.57 
o 

X 

X 

X 

X 

q 

HF 

free 

HF 

r =2.57 
o 

0 .0 

3.440 

3.274 


1.0 

1.844 

1.926 

0.1 

3.426 

3.261 


1.2 

1.473 

1.545 

0.2 ’ 

3 .380 

3.221 


1.4 

1.170 

1.221 

0.3 

3 .292 

3.148 


1.6 

0.931 

0.968 

0.4 

3.154 

3.045 


1.8 

0.754 

0.775 

to 

• 

O 

2.971 

2.906 


2.0 

0.616 

0.629 

0 .6 

2.756 

2.740 


3.0 

0. 276 

0.277 

0.7 

2.525 

2.550 


4,0 

0,160 

0.161 

0.8 

2.288 

2.342 


5.0 

0.106 

0.106 

0.9 

2.060 

2.134 
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Table 2.7 

2 

Compton Profile of O 


q 

HP 

free 

HF 

r =2.646 
o 

1 

1 q 

X 

HP 

free 

HP 

r =2.646 
o 

O 

♦ 

o 

5.426 

3.930 

o 

♦ 

% — i 

1.546 

1.870 

0.1 

5.358 

3.906 

1.2 

1.186 

1.414 

0.2 

4.892 

3.840 

1.4 

0.916 

1.06 2 

0.3 

4.140 

3.716 

1.6 

0.710 

0.814 

0.4 

3 .478 

3.530 

1.8 

0.566 

0.632 

0.5 

2.986 

3 . 286 

2.0 

0.464 

0.502 

0 .6 

2.596 

3.006 

3.0 

0.222 

0.226 

0.7 

2.280 

2.712 

4.0 

0.134 

0.140 

o 

• 

00 

2.004 

2.410 

■ 5.0 

0.088 

0.086 

0.9 

1.776 

2.126 
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Table 2>8 

Compton Profile of Cl 


q 

HP 

free 

HFS 

free 

0=1 

HP 

r =3.42 
o 

HPS 

r^=3 .42 
o 

o =1 .1 

HFS 

r =3.42 
o 

a =1 ,0 

HFS 

r =3.42 
o 

a =0.7 2 

0,0 

6.018 

5 .833 

5.796 

5.482 

5.568 

5.929 

0.1 

5.975 

5 .798 

5.760 

5.452 

5.531 

5.891 

0.2 

5.844 

5.664 

5.640 

5.355 

5.425 

5.757 

0.3 

5.586 

5.399 

5.410 

5.177 

5 .246 

5.494 

0,4 

5.193 

5.015 

5.080 

4.913 

4.961 

5.110 

0.5 

4.707 

4.557 

4.669 

4.573 

4.600 

4.633 

0.6 

4.176 

4.073 

4.204 

4.176 

4.182 

4.159 

0 .7 

3.654 

3.604 

3.736 

3 .756 

3.758 

3.668 

0.8 

3.178 

3.173 

3.292 

3.335 

3 .321 

3.202 

0.9 

2.748 

2.790 

2.882 

2.944 

2.929 

2.776 

1 .0 

2.405 

2.460 

2.523 

2.592 

2.571 

2.108 

1.2 

1.876 

1.946 

1.957 

2.028 

2.014 

1.865 

1.4 

1.530 

1.593 

1.578 

1.637 

1.6 23 

1.5 29 

1.6 

1.306 

1.353 

1.331 

1.376 

1.370 

1.307 

1.8 

1.166 

1.188 

1.172 

1.201 

1.202 

1.154 

2.0 

1.051 

1.070 

1.059 

1.073 

1.07 2 

1.054 

3.0 

0 .7 28 

0.732 

0.7 27 

0.732 

0.729 

0.730 

4.0 

0.492 

0.504 

0.496 

0.504 

0.502 

0.502 

5.0 

0.3 25 

0.364 

0.327 

0.365 

0,364 

0.364 
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Table 2.9 



Compton 

Profile of 

Br~ 




r~ 




HFS 

i 


HFS 

q 

free 

i 

q 

free 



1 


a=l 

0.0 

7.998 


1*0 

3.946 

0.1 

7.961 


1.2 

3.406 

0.2 

7.711 


1.4 

3.054 

0 .3 

7.367 


1.6 

2.816 

0.4 

6.959 


1.8 

2.645 

0.5 

6.395 


2.0 

• 2.506 

0.6 

5.793 


3.0 

1 . 884 

0.7 

i . 

5.229 


4.0 

1.328 

0,8 

4.730 


5.0 

0.944 

0,9 

4.306 
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shown in the Tables 2.10 through 2.20. Also shown in these 

tables are the results of Aikala^^ using the overlap model. In 

« 

all the Cases examined/ we see that the J{0) value is overesti~ 
mated by both the HP and HFS («=!) free-ion CP's as compared to 
the corresponding experimental CP's. Thus the model employing 
free-ion CP' s is not a good one to study the CP' s of ionic crys- 
tals. It may also be noted here that the agreement between HP 
and HFS (“=l) inner-shell CP' s of free ions is excellent in all 
the cases examined. However, the valence-CP' s in the HP and HFS 
(o=l) frameworks differ considerably. 

We will now examine the HF-crystal-ion CP' s presented in 
Tables 2.1 through 2.9. We see that in the case of all the 
crystal-anions examined, viz., F and 0 and Cl , the J(0) value 

is less than that in the corresponding free HF anions. On the 

+ 24* 4" 

other hand, in the case of crystal-cations, viz., Li , Be , Na , 
Mg and K , the HF crystal-ion J(0) value is more than the corres- 
ponding free ion J(0) value. The external spherical-potential 

13 

used to simulate the behaviour of ions in crystals is known to 
expand the charge distribution in the case of cations and con- 
tract it in the case of anions. From this behaviour, it may be 
expected that in the case of anions, the crystal ion CP* s should 
be broader than the corresponding free-ion ones, leading to a 
lower J(0) value than that for the free ion. On the other hand, 
in the case of cations, the crystal-ion CP* s are expected to be 
narrower than the corresponding free ion CP. These expected 
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Table 2.10 



Theo retical 

and Experimental 

CP' s of LiF 


q 

HF 

free 

HF 

crystal 

Overl ap 
model 

Expt. 

0.0 

4.096 

3.948 

3.855 

3.849 

0.1 

4.078 

3.933 

3.840 

3.822 

0.2 

4.024 

3 .831 

3.794 

3.776 

0.3 

3.920 

3.792 

3.716 

3.704 

0.4 

3 .762 

3,665 

3.603 

3.608 

if ) 

• 

o 

3 .555 

3 .550 

3.461 

3.478 

0.6 

3 .312 

3.302 

3.293 

3.318 

0.7 

3.052 

3.078 

3.101 

3.124 

o 

00 

2.782 

2.838 

2,886 

2.901 

0.9 

2.5 20 

2.592 

2.653 

2.666 

1.0 

2.270 

2.348 

2.411 

2.427 

1.2 

1.833 

1.899 

1.932 

1.951 

1.4 

1.468 

1.513 

1.517 

1.519 

1.6 

1.175 

1.206 

1.197 

1.178 

1 .8 

0.952 

0.969 

0.96 2 

0.942 

2.0 

0.776 

0.785 

0.787 

0.769 

3.0 

0.330 

0.329 

0.330 

0.341 

4.0 

0.178 

0.179 

0.181 

0.19 2 

5 .0 

0.116 

0,116 

0,113 

0.112 
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Table 2.11 



Theoretical 

and Experimental 

CP' s of NaF 


q 

HF 

free 

HF 

Crystal 

Overlap 

model 

Expt. 

0.0 

5,725 

5.568 

5.559 

5.662 

0.1 

5.708 

5.549 

5.542 

5.622 

0.2 

5.647 

5.495 

5.488 

5.531 

0.3 

5.532 

5.394 

5.391 

5.395 

0,4 

5.364 

5.261 

5.250 

5.231 

0.5 

5.089 

5.030 

5.068 

5.044 

0.6 

4.860 

4.850 

4.850 

4.828 

0.7 

4.565 

4.590 

4,599 

4,586 

0.8 

4.250 

4.303 

4.322 

4.324 

0.9 

3.942 

4.016 

4.026 

4.045 

1.0 

3.636 

3 .718 

3.723 

3.753 

1.2 

3.075 

3,147 

3.130 

3.175 

1.4 

2.574 

2.625 

2.605 

2.6 23 

1.6 

2.149 

2.180 

2.167 

2.132 

1.8 

1.802 

1 .617 

1.808 

. 1.750 

2.0 

1.512 

1.519 

1.516 

1.473 

3 .0 

0.700 

0.699 

0.706 

0.700 

4 .0 

0.392 

0.393 

0.396 

0,401 

5 .0 

0.254 

0.254 

0.253 

0. 241 
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Table 2.12 


Theoretical and Experimental CP* s of KP 


q 

HF 

free 

HF 

crystal 

Expt. 

0,0 

7.868 

7.730 

7.654 

0.1 

7.832 

7.699 

7.619 

0.2 

7.738 

7.608 

7.507 

0.3 

7.5 64 

7.443 

7.321 

0.4 

7.297 

7.213 

7.086 

0.5 

6.949- 

6.902 

6.809 

0.6 

6 .530 

6.529 

6.489 

0.7 

6.066 

6.096 

6.106 

0.8 

5.576 

5.6 24 

5.694 

0.9 

5.084 

5.15 2 

5.265 

1.0 

4.606 

4.674 

4.785 

1,2 

3.748 

3.796 

3 . 890 

1 .4 

3.034 

3.079 

3.127 

1.6 

2.481 

2.508 

2.535 

1.8 

2.070 

2.086 

2.095 

2.0 

1.762 

1.774 

1.780 

3 .0 

1.023 

1.023 

1.021 

4.0 

0.701 

0.702 

0.683 

5.0 

0.497 

0.496 

0.487 
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Table 2.13 


Theoretical and Experimental CP‘ s of LiCl 


q 

HF 

free 

HFS 

free 

o(=l 

HF 

crys tal 

HFS 

crystal 
a =1 

Overlap 

model 

Expt. 

0,0 

6.674 

6.497 

6.470 

6.284 

6 . 245 

6 .360 

0.1 

6.6 27 

6.458 

6.432 

6.243 

6.209 

6 .207 

0.2 

6.488 

6.314 

6.300 

6.125 

6.096 

6.163 

0.3 

6.214 

6.035 

6.054 

5.926 

5.904 

5 .932 

0.4 

5.801 

5.629 

5.700 

5.618 

5.642 

5.6 25 

0.5 

5.291 

5.137 

5.263 

5.224 

5.312 

5.264 

0.6 

4.732 

4.633 

4.766 

4.772 

4.911 

4.861 

0.7 

4.180 

4.132 

4.264 

4.310 

4.443 

4.421 

0.8 

3.672 

3.669 

3.786 

3.833 

3.935 

3.962 

0.9 

3.208 

3.250 

3.340 

3.401 

3.432 

3.506 

1.0 

2.831 

2.886 

2.945 

3.003 

2.973 

3.07 2 

1.2 

2.236 

2.305 

2.311 

2.268 

2.270 

2.343 

1.4 

1.828 

1.889 

1.870 

1.909 

1.836 

1.854 

1.6 

1.550 

1.595 

1.619 

1.598 

1.562 

1.550 

1.8 

1.364 

1.384 

1.366 

1.382 

1.366 

1.357 

2.0 

1. 211 

1.226 

1.215 

1.214 

1 . 217 

1.221 

3.0 

0.782 

0.784 

0.779 

0.773 

0.780 

0.775 

4.0 

0.510 

0.524 

0.514 

0.516 

0.513 

0.493 

5.0 

0.333 

0.372 

0.335 

0.370 

0.337 

0.334 
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Table 2.14 

Theoretical and Experimental CP* s of NaCl 


q 

HP 

free 

HFS 

free 

HP 

crystal 

HFS 

crystal 
a =1 

Overlap 

model 

Expt. 

0 .0 

8.303 

8.069 

8.090 

7.896 

8.001 

8.130 

0.1 

8.257 

8.034 

8.048 

7.855 

7.961 

8.040 

0.2 

8.112 

7.886 

7-914 

7.735 

7.834 

7.861 

0.3 

7.826 

7 .595 

7 .656 

7.534 

7.611 

7.601 

0.4 

7.403 

7.178 

7 . 296 

7.211 

7.294 

. 7.263 

0.5 

6.825 

6.693 

6.793 

6.801 

6.893 

6.858 

0.6 

6.280 

6.136 

6.314 

6.317 

6.412 

6.404 

0.7 

5.694 

5.603 

5.776 

5.823 

5.871 

5.919 

0.8 

5.140 

5.098 

5 .253 

5.302 

5.307 

5.417 

0.9 

4.630 

4.632 

4.764 

4.813 

4.762 

4.921 

1 .0 

4 . 197 

4.215 

4.315 

4.358 

4.271 

4.455 

1.2 

3.478 

3.518 

3.559 

3.589 

3.494 

3.630 

1.4 

2.934 

2.978 

2.982 

2.990 

2.944 

2.979 

'1,6 

2.524 

2,561 

2.543 

2.548 

2.531 

2.493 

1.8 

2.204 

2.236 

2.214 

2.220 

2.206 

2.166 

2.0 

1.957 

1.971 

1.949 

1.948 

1.948 

1.910 

3.0 

1.152 

1.175 

1.149 

1.148 

1.149 

1 . 160 

4.0 

0.724 

0.750 

0,7 28 

0.738 

0.726 

0.706 

5.0 

0,471 

0.516 

0.473 

0.514 

0.476 

0.453 
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Table 2.16 



Theoretical 

and Experimental CP' s of MgO 


q 

HF 

free 

HF 

crystal 

Overlap 

model 

Expt. 

0.0 

7.400 

5.932 

5.597 

5.840 

0.1 

7.3 28 

5.906 

5.575 

5-821 

0 . 2 

6.852 

5.830 

5.507 

5.758 

0.3 

6.094 

5 .692 

5.39 

5.641 

0.4 

5.402 

5.472 

5.270 

5.465 

0.5 

4.870 

5.204 

5.127 

5 . 23 2 

0.6 

4.454 

4.888 

4.946 

4.951 

0.7 

4.096 

4.552 

4.705 

4.632 

0.8 

3.774 

4.198 

4 . 400 

4. 290 

0.9 

3 .490 

3 .942 

4.048 

3.938 

1.0 

3.196 

3.538 

3.678 

3.595 

1.2 

2.708 

2.948 

2.891 

2.966 

1.4 

2.300 

2.446 

2.432 

2.425 

1.6 

1.95 2 

2.050 

2.021 

1.988 

1.8 

1.670 

1.724 

1.712 

1.665 

2.0 

1.436 

1.46 2 

1.467 

1.414 

3.0 

0.500 

0.492 

0.739 

0.710 

4.0 

0.304 

0.310 

0.411 

0.409 

5.0 

0.204 

0 . 208 

0.259 

0.250 
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Table 2.17 

Theoretical and Experimental CP' s of BeO 


q 

HP 

free 

HP 

crystal 

Expt, 

0.0 

5.898 

4.450 

3.80 

0.1 

5.830 

4.424 


0.2 

5.360 

4.354 

3.73 

0.3 

4.610 

4.220 


0.4 

3.934 

4.022 

3.54 

0.5 

3.432 

3.764 

- 

0.6 

3.030 

3.468 

3.22 

0.7 

2.702 

3.156 

- 

0.8 

2.410 

2.834 

2.76 

0.9 

2.168 

2.528 

- 

1 .0 

1.922 

2.250 

2.27 

1,2 

1.528 

1.752 

1.81 

1.4 

1.222 

1.358 

1.45 

1 .6 

0.982 

1.07 2 

1.15 

1 .8 

0.804 

0.854 

0.93 

2.0 

0.67 2 

0.694 

0.79 

3 .0 

0.318 

0.314 

0.43 

4 .0 

0.178 

0.176 

0.25 

5 .0 

0.108 

0.106 

0.14 
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Table 2.18 

Theoretical and Experimental CP' s of LiBr* 


q 

HFS 

free 

ct=l 

Overl ap 
rrodel 

Expt. 

0.0 

8.59 2 

8.320 

8 . 265 

0,1 

8.551 

8.274 

8-191 

0.2 

8.291 

8.129 

8,000 

0.3 

7.933 

7.886 

7.705 

0.4 

7.502 

7.557 

7.327 

0.5 

6.905 

7.142 

6.890 

0.6 

6 . 283 

6.630 

6.422 

0.7 

5.687 

6.037 

5.939 

0.8 

5.156 

5.415 

5.449 

0.9 

4.697 

4.830 

4.974 

1.0 

4.398 

4.331 

4.529 

1,2 

3.697 

3.638 

3.789 

1.4 

3.282 

3.239 

3.288 

1 .6 

2.990 

2.974 

2.963 

1.8 

2.774 

2.767 

2.7 26 

2.0 

2.596 

2.592 

2.551 

3.0 

1.870 

1.850 

d.889 

4.0 

1 . 280 

1.254 

1. 266 

5.0 

0.884 

0.854 

0.843 


*ls Gontribution to CP of Br~ subtracted out in order to 
enable comparison with literature data ( ref . lo) . 
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Table 2.19 



Theoretical and Experimental CP' s of 

NaBr* 


q 

HFS 

free 

Overl ap 
model 


Expt . 

0 .0 

10.164 

10.093 


10.053 

0.1 

10.127 

10.043 


9.969 

0.2 

9.863 

9.882 


9.779 

0 .3 

9.494 

9.599 


9.483 

0 .4 

9.05 2 

9.205 


9.078 

0.5 

8.443 

8.705 


8.585 

0 .6 

7.786 

8.106 


8.027 

0.7 

7.158 

7.442 


7.436 

0.8 

6.585 

6.774 


6.839 

0.9 

6.078 

6.161 


6.270 

1 .0 

5.631 

5.640 


5.76 2 

1.2 

4.978 

4.873 


4.938 

1.4 

4.369 

4.346 


4.334 

1 .6 

3.956 

3.941 


3.897 

1.8 

3.628 

3.609 


3.586 

2.0 

3.339 

3.395 


3.298 

3.0 

2.259 

2.218 


2.288 

4.0 

1.506 

1.467 


1.465 

5 .0 

1.030 

0.993 


1.005 

* Br~ Is core contribution 

comparison with literature 

subtracted out in 
data (ref. lO) , 

order 

to enable 



69 


Table 2,20 

Theoretical and Experimental CP* s of KBr* 


q 

HFS 

free 

Overlap 

model 

Expt. 

0.0 

12.405 

12.238 

12.347 

0.1 

12.350 

12.176 

12.265 

0.2 

12.057 

12.978 

12.030 

0.3 

11.636 

11.630 

11.620 

0.4 

11.103 

11.135 

11.086 

0.5 

10.360 

10.503 

10.461 

0 *6 

9.525 

9.752 

9.767 

0 .7 

8.701 

8.927 

9.015 

0 .8 

7.933 

8.097 

8.244 

0.9 

6.600 

7.322 

7.497 

1 .0 

5.962 

6.639 

6.805 

1.2 

5.619 

5.564 

5.683 

1 .4 

4.806 

4.806 

, 4.922 

1 .6 

4.271 

4.271 

4.351 

1.8 

3.886 

3.882 

3.893 

2,0 

3.579 

3.578 

3.563 

3,0 

2.572 

2.541 

2.611 

4.0 

1.807 

1.777 

1.748 

5.0 

1.264 

1.233 

1.151 


*Gontribution to CP from Is shell of Br stibtr acted out in 
order to enable comparison with literature data (ref. lo) . 
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trends are borne out by the HF free and crystal“ion CP data for 

2-f 4- 2+ + ~ 2— *- 

Li , Be / Na , Mg , K , r , 0 and Cl . We may note here 
that the CP's of inner shells are nearly identical for each of 
the HF free and crystal-ions examined. The CP' s of ionic crys- 
tals LiF, NaF, KF^ LiCl, NaCl, KCl, BeO and MgO (sho-wn in 
Tables 2.10, 2.11, 2.12, 2.13, 2.14, 2.15, 2.19 and 2.20 res- 
pectively) were calculated from the HF crystal-ioh data. These 
crystal CP' s can be seen to be in better agreement with the 
experimental ones than the corresponding f ree-ion CP' s. These 
Tables also show the CP's calculated for these ionic crystals 
by Aikala et al.^^ with a more elaborate overl ap model . The 
CP's calculated by the overlap model give, in general/ better 
agreement with the experimental ones. 

We shall now discuss the main features shown by HFS crystal- 
ion CP's. Firstly we note that the inner-shell crystal-ion CP's 
do not differ significantly from the corresponding free-ion onfes. 
Secondly, J(0) value for crystal-cation is more than that for the 
corresponding free cation, whereas J(o) value for crystal-anion 
is less than that for the free anion. This is in agreement with 
the earlier discussion of Paschal is and Weiss in terras of contrac- 
tion of anions and expansion of cations when placed in an exter- 
nal spherical potential. Thirdly/ the HFS crystal-ion CP for a 
given ion varies significantly with a . This is seen to be the 
case for all the ions examined (viz., Li , Be , Na , Mg , K 
and Cl”) . We notice that in all the ions under study, the J(0) 
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value decreases with increasing a , This behaviour can also be 

expected since a numerically high value of exchange potential 
. 19 

is known to lead to a contracted charge distribution* S.uch a 
contracted charge distribution implies a broader CP and lower 
J(0) value. Examining the effect of variation of a on the 
CP's orbital wise reveals an interesting phenomenon. The 
valence orbital CP's are more saasitive to a -variation than the 
core orbital CP's. To illustrate this we list the different 
orbital CP* s in Mg^"'' in the HFS framework with the Watson sphere 
potential as a function of three different a values (see Table 2.2l)\ 

Considering a purely as a parameter in the spirit of Lindgren 
5 6 

and Schwarz * we see from the above tables^ that there should 

be a good agreement between the HF and HFS crystal-cation CP' s 

when a value of a > 1 is used in the HFS calculations. The J(C) 

“i" / 

value for K calculated from HFS crystal don wavefunction with 
“=1.1 can be seen to be in good agreement with the corresponding 
HF crystal ion J(o) value. However# in the case of Cl ion# as 
seen from Table 2,8# a value of ot lying between 0.72 and 1.0 is 
expected to give good agreement between the HFS and HF crystal- 
ion CP's. It is clear from these results that with the addition 
of the external Watson sphere potential to the ionic HFS hamil- 
tonian the nature of exchange-correlation potential is con side- 
rably altered from that of the free ion. This results in the 
necessity for correcting the exchange correlation potential in a 
manner which is different for the core and valence orbitals. 
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The CP's of ionic crystals calculated from the HFS (a=l) 
crystal-ion CP's are shown in Tables 2.12, 2.13 and 2.14 for 
LiCl, NaCl and KCl respectively. . These CP' s are, in general, in 
better agreement with the corresponding experimental CP' s than 
the respective free-ion HFS (a=i) CP's. However, HFS crystal 
ion results with a =1 do not appear to be as good as HF crystal 
ion results. Better agreement might result in HFS calculations 
vAen a is treated as a parameter. 

Commenting on the work of Paschalis and Wei ss^^ Weyrich^^ 
has pointed out that the virial theorem is violated in their model 
employing external crystalline potential due to two reasons: 

(i) The constant potential inside the sphere does not 

obey — 2 force-law. 
r 

(ii) The STO' s used in their HF calculation were not opti- 
mized in the case of crystal-ion wavefunction. 

We may note here that the Watson model using HFS wavefunc- 
tions does not suffer from the defect (ii) above due to the 
numerical procedure employed in obtaining the wavefunctions. 

This fact, coupled with the approach of using a as a paraneter 
makes the HFS-crystal ion model attractive. 

II. 4 SUMMARY AND CONCLUSIONS 

CP' s were calculated from HF and HFS (a=i) free ion wave- 
functions as well as HF and HFS (with “-variation) crystal-ion 
wavefunctions employing an external spherical potential proposed 
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by Watson. Use of these crystal-ion CP's gives better agreement 
with experiment than the use of corresponding free ion CP' s- This 
result holds good both in the HP as well as HFS frameworks. The 
HFS crystal-ion CP's for valence orbitals show a marked depen- 
dence on the value of the a used. Treating a purely as a para- 
meter^ it has been shown that a value of a greater than unity 
gives better agreement between HP and HFS crystal-ion CP' s for 

-j- 0“i" “T* 0“h ”4“ 

all the cations examined (viz., Li , Be , Na , Mg and K ) 
while, for Cl , such an agreement may be expected for a value 
of a less than unity. The introduction of an external potential 
in the free ion HFS hamiltonian alters the delicate representation 
of the exchange-correlation term. 

One of the reasons for the violation of virial theorem in 
the case of HP crystal-ion wavefunctions is that the STO exponents 
used in these wavefunctions have not been optimized. Due to their 
numerical nature HFS crystal ion wavefunctions do not suffer from 
this defect. HFS crystal ion wavefunctions obtained with para- 
metrized a seem to offer promise in the calculation of crystal 


CP' s 
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CPiAPTER III 

ELECTRON MOMENTUM DISTRIBUTIONS 
ERDM valence-bond WAVEFUNCTIONS* 


III.l INTRODUCTION 

In 1941 Coulson^ presented the first paper on the chemi- 
cal interpretation of electron momentTm distributions (EMD) in 
molecules. He analyzed the EMD' s obtained from the Heitler- 
London (HL) as well as molecular orbital (MO) wavefunctions for 

diatomic molecules. The momentum- space (p-space) wavef^anctions 

2 

were obtained by him by Dirac transformation of the correspond- 
ing co-ordinate-space (r-space) wavefunctions; 


X ( * ^2 ^ * ^N^ ^ ^ 


-m f 

w)= J 


-i{ Pt . r, + p„ . r^ t + p„- r„ } 

5 • Y v 2 2 




. .. (3.1) 



* The material presented in this chapter has been accepted for 
publication in the form of a paper; S.R. Cadre and P.T, 
Narasimhan# Int. J. Quant. Chem. (to be published) . 



where the electrons are numbered 1,2...N. Coulson also pointed 
out that the form of the wavef unction is retained in the Dirac 
transf ormation • ( 3 . 1 ) . Thus if one is working within the orbital 
approximation and with the wavefunction 




■ * 


(3,2) 


then it is easy to see on substituting in (3*l) that 


= X,(^) X2(p,)... X^(p^) 


(3,3) 


Similarly for the electron-pair method, where 




* ‘ ’^2N-1, 2N^^2N-1' ^2N^ 

' ^ V-\ 


(3.4) 


one obtains, on substitution in (3.1), 


^ ' *£2N^ “^1,2^^'P2^* ^3,4*‘^"^4^ ^2N-1, 2N^P2N-1' ^2N^ 


(3.5) 


The momentum densities for the HL and MO wavef\anctions v/ere 
obtained by Coulson by squaring the p-space wavefunction. Thus, 
for a bond A-B, the Hei tier- London wavefunction can be written as 
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where c})^ and <|)^ are the atomic orbitals (AO) centred at the 
nuclear positions and ^^ respectively and S is the overlap 
integral: 


Ur) 

a 


.U(r) 

jd 


dr 


(3.7) 


The following expression was obtained for a homonuclear 
diatomic molecule starting from the corresponding HL wavefunc— 
tion (Eq. 3.6) on Fourier transformation and subsequent squaring: 


l^(Pi/P2^ I' 


1 + 


cos{ (p, -p„).(r -r, )} 
c^2 ■.>^a Ob 


1 + S‘ 


A(Pj^) 


A(p,) 


(3.8) 


where p- and p- are the momenta of electrons 1 and 2 respectively 
and A(p ) is the Fourier transform of the AO (r) , viz., 

A(p) = (> (r) dr (3.9) 

Coulson pointed out that the momentum density corresponding to 
electron 1 can be obtained from (3.8) by integration over P 2 . 

It can also be obtained more easily by employing a mixed wave- 
function of p. and r^ followed by integration over the co-ordi- 
nate r^- The momentum densities obtained from valence-bond (VB) 
wavefunctions by these two procedures can be shown to lead to 
identical results. This technique of using a mixed wavef-unction 
to obtain momentum density corresponding to electron 1 was 



1 

Introduced by Coulson who obtained the following expression 
for this quantity: 


X ( ) 1 


1 + S cos { p. . ( r -r, ) > 

V/-»3 v-r«lD 


lA(p^) 


(3.10) 


The corresponding result for a MO wavefunction for a homo- 
nuclear diatomic molecule 


9 1/2 

(r) = ( *? t (>,)/( 2+- 2S ) (3.11) 

%rv 3D 

is 

|x (p) P 

It was pointed out by Coulson that the momentum densities in 
(3.10) and (3.12) show maxima when the electron moves perpendi- 
cular to the bond direction. Expression (3.8) shows additional 
maxima for p, = p,p. This means, according to Coulson, that if 
atoms A and B are considered in terms of oscillators/ there is 
binding" when the oscillators are in phase • Considering the 
repulsive-state wavefunction 


1+ cos{p. (r -r, ) } 


(3.12) 


he obtained 
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and by employing mixed wavef unction involving and r„ followed 
by integration over r^, the niomentiirn density corresponding to 
electron 1 is obtained as 




1-S cos { Pt . (r -r, )} 

V^a v"vL) 

1 - 


|a (^) 


(3.15) 


According to (3.14)/ the momentum density vanishes .for 



Expression (3.15) shows a minimum when the electron moves per- 
pendicular to the bond. 


In this chapter/ we present an analysis carried out by us 
on lines similar to those of Coulson but explicitly including 
ionic terms as well. To the best of our knowledge/ no such 
analysis has been done previously. Earlier, Hicks^ calculated 

the shape of the Compton line for using VB wavefunction of 

4 . 

Weinbaum which includes ionic terms as well as polarization 

(2p orbital). However, Hicks did, not analyse the detailed 
characteristics of VB momentum densities and their dependence on 
ionic terms. Here we will discuss specifically the effect of 
adding ionic terms on the momentum densities (Section III. 2). In 
Section III. 3 we discuss the Compton profile (CP) of H 2 with a 
view to assess the importance of ionic terms as well as polari- 
zation. We also present the CP calculated from an MO wavefunction 
of Shull and Ebbing^ which employs floating functions . Comparison 
with other avail^l© theoretical data as well as experimental CP 



83 


is also made. In Section II 1, 4, the VB wavef unction constructed 
from the orthogonalized atomic orbitals (OAO) is examined in the 
p—space. The results are then compared with an analysis of 
Slater^ who showed that an HL wavefunction constructed from OAO' s 
shows no minimum diagonal energy for the H 2 molecule. 


III. 2 WAVEFUNCTION WITH IONIC- TERMS ^ 

The r-space wavefunction for a homonuclear diatomic mole- 

— .ill. 

cule A-B including the ionic teinns A B and A B is 


’J' (r. , r^) 


<:>^(1)<>, ( 2) +0, (1) (2)+ x{<) (!)<' {2)+ <>, (!)<>, ( 2) } 

D D a a a d d 

{ ( 1 + X ^) ( 2+2S^) + 8X S} 


(3.16) 


This yields the following expression for the momentum density 

|A(pJ 1 lA(p ) 1^ 2 

I x(Pt,P,) P = — ^ l-*-co/(Pi-P2^ • (^="5^33^ ^ 

' ^12' (1+x2)(i + s2)+4Xs 

+ X^cosi(£i+P2).(r^-5,) 3 


+ X (cos[Pj . (r-rj,)] + cos [ Pjt Pa-ib’]” 

(3.17) 


On comparing this result with Eq. (3.8) for the purely 
covalent HL wavefunction (3.6), we find additional terms such 



factor". Hence no definite phase relationship which leads to 
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2 

maxima in j x^P^PoM can be obtained in terms of the oscilla- 
tor picture. Thus the "in-phase" oscillator picture seems to 
be valid only for the purely covalent HL wavef unction. We also 
note here that the purely ionic terms A^B and A individually 
give rise to the term cos{(p +p„). (r -r, ) } which implies that 
the momentum density is maximum for p^= - P 2 « Thus the momentum 
density for a purely ionic term shows a maximum when the oscilla- 
tors are "out-of-phase". The momentum density corresponding to 
electron 1 obtained from the VB wavefunction (3.16) following 


the use of a mixed wavefunction involving p^ and r 2 is: 


X (pJ 


l+X^t 2XSt {S+SX^+ 2X}cos{P-, . (r -r, )} „ 

5 ; lACpUl^ 

{(1+X^) (1+S^)+ 4 XS} 


(3,18) 


Analogous to (3.5) and (3.8), this expression also shows maxima 
for p, .(r “r, ) =0. Thus the present results obtained for VB 
wavefunction (3.11) agree with the earlier analysis by Coulson, 
of HL and MO wavefunctions, in terms of greater probability of 
an electron moving perpendicular to the bond rather than along 
it. However, the phase relationship, in the oscillator picture, 
existing for the HL wavefxanction is seen to be aibsent for the VB 
wavefunction which includes "covalent as well as ionic terms. 
It must, however, be borne in mind that this phase relationship 
is an artefact of the purely covalent VB wavefunction and should 
not be considered to have any physical significance. 
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We will now present a similar analysis for the waveftmc- 
tion of a heteronuclear diatomic molecule. . If# according to 
chemical intuition based upon the electronegativities/ only one 
ionic teim, say A B , is included then the expression for the 
momentum density apart from the normalisation constant is, with 
an obvious notation/ 


1 


x(p^ , p^) 1 


^ = 1 A(p ) PIbCp + 1a(pJ12. jB(p;) l2 + 

x2lA(p^)j2, 1a(p2)| ^ + 2 jA(^)j.l B(p^)j . Ia (P2) I . jB<^)i . 



(3.19) 

If we take into account the fact, as previously pointed out by 

1 ■ - ^ 

Coulson that in most cases where bonds are formed, conditions 

of symmetry imply arg A{p) = arg B(p), the above expression is 
considerably simplified. The result thus obtained from (3,19) 
by imposing this symmetry criterion is similar to (3.17), but the 
term cos {(p-+p.,).(r -r,) } is absent, corresponding to the inclu- 
sion of only one ionic term A~B , Such a term reappears on inclu- 
Sion of both the ionic terms A^ B and A B . The general conclu- 
sion that can be drawn from (3*19) is similar to that obtained 
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from (3. 17 )/ viz., the "phase relationship" in the oscillator 
picture becomes complex with the introduction of one or more 
ionic terms even in the heteronuclear cases. 


III. 3 EFFECT OF INCLUSION OF IONIC TERMS MJD POLARIZATION 
ON THE COMPTON PROFILE OF 

3 

As mentioned earlier in Section III.l, Hicks calcu- 

4 

lated the shape of the Compton line of using Weinbaum' s 
wavefunction constructed from AO's of the form ls(r)+ cT2p(r) 
and including ionic terms as well ( wave f -unction of ionic-polari- 

n 

zation type). More recently Eisenberger has reported the CP 
of H 2 measured with 17.37 keV Mo X-rays, alongwith the theore- 
tical HF-SCF and MC-SCF CP's. We may note here that the agree- 
ment between theoretical and experimental CP* s is not good even 

s 

with high quality MC-SCF wavefimction. Brown and Smith"* have 
calculated CP of H 2 using Liu* s configuration-interaction (Cl) 
wavefunction which includes 39 configurations. This gives 
satisfactoiry agreement with experiment. We may note here that 
the comparison between theory and experiment must be done cauti- 
ously due to several experimental factors. The experimental 

■9 ■" 

results have to be corrected for multiple scattering. However, 
multiple scattering is not expected to be appreciable in the 
c ase of H 2 • The deconvolution procedure employed to remove the 
effects of instrumental broadening is known to lead to oscilla- 
tory behaviour^ ^ and yet the results have to be corrected for 
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the residual instrumental function. 

Since the ionic terms in VB theory play a role similar to 

that of Cl in MO theory, it is desirable to analyze the effect 

of adding ionic terms to the HL wavefunction on the CP of H 2 . 

1 1 

We have examined here the simple VB wavefunctions of Wang, 

4 S 

Weinbaum and MO wavefunction of Shull and Ebbing employing 

floating functions. The details of these wavef\inctions are 

3 

given in Table 3,1. Hicks calculated the shape of the Compton 

4 

line using Weinbaum s wavefunction of ionic-polarization type 
but assuming experimental intemuclear distance of 1.40 a.u. He 
carried out the integrations involved in the calculation graphi- 
cally to within 1 % accuracy. We have repeated the calculations 
of Hicks for the intemuclear distances of 1,40 a.u. and 1.43 a.u. 
The momentiom densities & CP' s given by these two sets of values 

agree with each other to within 0.2%. The CP's, within the 

12 

impulse approximation, are obtained from one-electron momen- 
tum densities via (3.20) and (3.21). 


I(p) = Jp^P(p) sin ep. depd<>p 


(3.20) 


J(q).= j 

1^1 


I(p) 

2p 



(3.21) 


For the wavefunctions employing only s-type AO's, the 


integration in( 3. 20) ,is simple. Though more 'complex,. , it can be. done 
analytically for a wavefunction involving polarization. The 



Details of the Wavefunctions Employed in the Present Compton Profile Calculations 


88 





89 


integration in (3*21) was carried out numerically by using 
Simpson s rule. The CP* s obtained were normalized to two 
electrons by 

J(q) dq = 1 (3.22) 

O'' 

Overlap integrals were evaluated by the standard formulae of 
1 3 

Mulliken et al. Theoretical CP, j(q), normalized to two elec- 
trons for the above mentioned wavefunctions alongwith the experi- 
mental and Cl-thsoreticai profiles are presented in Table 3,2. 

The CP' s calculated from Weinbaxm' s wave fxmct ions can be seen 
to be quite close to the Cl-theoretical profiles. The J(0) value 
for the polarized wavefunction is higher than that for the unpo- 
larized one since the 2p orbital leads to a narrower CP. On the 
other hand/ the CP computed from Wang's wavefunction does not 
agree so well with either the experimental or Cl-theoretical 
profile in the smaller q- region. The observation that the CP 
calculated from VB wavefunction including ionic terms agrees 
well with the Cl-theoretical profile in the case under study 
does not imply that the CP's obtained from the VB wavefunctions 
will always give a very good agreement with the corresponding 
Cl— theoretical profiles. However/ the addition of ionic terms 
seems to improve the EMD' s markedly. In the light of this 
experience with H 2 molecule one may say that the inclusion of 
ionic terms would be quite necessary for the calculation of CP s 
of heteronuclear molecules from VB wavefunctions. 
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1 1 1. 4 ANALYSIS OF OAO RESULTS O F SLATER IN TERMS OF MOMEt-ITUM 
DENSITIES ~ 

6 

Slater has shown/ with the example of molecule that if 
a VB calculation is done using only the covalent wavef\inction 
constructed from OAO' s instead of overlapping AO' s, the resultant 
diagonal energy of the lowest singlet shows no minimum. We present 
here our momentum density results from a wavefunction constructed 
from OAO' s and ^ which are/ in turn/ linear combinations of 
the overlapping AO's and'I’j^ 






( 3 . 23 ) 


where and C 2 are determined by the normalization condition 


Ci^ + 2C^C2S + = 1 


(3.24) 


S is the overlap between and<)j^ and the orthogonality 
between and is : given by 


f Oa Ob ° 

j 

which leads to . - 

2 C,Cj+ (C^^+Cj^ls >0 


( 3 . 25 ) 


( 3 . 26 ) 



93 


Following Slater, 
the OAO' s is 


the HL wavef unction constructed from 


’I'cov = t «^(2) ^^(1) 


1 

V2(l-S^) 


<^(1)<I>^(2) >1 


S{« (1) f) (2) + 

a 

(3.27) 


With this wavefunction, the momentum density corresponding to 
electron 1 is 



1-S cos {Pt( 




) } 


1-S' 


1 A(p^) 


(3.28) 


It will be seen that this equation is exactly identical 
with Eq. (3.15) obtained for the repulsive-state wavefunction 
(3.13) . Thus, since the momentum densities corresponding to 
electron 1 are identical, for all the values of p^ , in the case 
of wavefunctions (3.27) and (3.13) and as (3.13) is a repulsive 
state, (3.27) also may not be expected to show binding. This 
conclusion is in accordance with the earlier discussion of Slater. 

We can reach the same conclusion on comparing the "diffrac- 
tion factors" for (3.6) and (3.27) . The diffraction factor for 
the former is maximum whereas that for the latter is minimum when 
the electron moves perpendicular to the bond. The latter reaches 
maximum for p^.(r -r. ) -r where jA(p 2 )| is already quite small. 
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Thus, for the wavefunction (3.27), the probability of the elec- 
tron moving perpendicular to the bond is less than that for the 
motion along the bond direction and hence it should not show 
binding. 


The ionic wavefunction constructed from OAO' s is 


^onic= 


(3.29) 


f2(l-S^) 


(2)+ « (2) 4 (1) O 


This wavefunction leads to an expression for the momentum 
density which is identical to (3.28). Thus, purely covalent and 
ionic wavef unctions constructed from 0A0‘ s yield momentum densi- 
ties identical to that for the repulsive-state wavefunction. 

Wavefunctions (3.13), (3.27) and (3.29) can be shown to lead to 

14 

identical r-space density matrices of first order. 


Y (^ I r^^) where 



’J' (r-, , r„) 


d r, 

V.A,* 


(3.30) 


The corresponding density matrix in the iromentum space is 
given by 


Y (P' l P) 

Sdtfs I 


= (2Tr)”^/ 


(r' 


Ir) 


-i(p.r-p' .rV) 


V% ^ 


dr dr’ 


(3.31) 
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Since Y'( r' 1 r) is identical for these wavefinctions/ .Y(p' 1 p) 

is also identical. As discussed above/ for a given value of the 

momentum/ the momentum density perpendicular to the bond direc~ 

tion is less than that along the bond direction and hence these 

wavef unctions are not expected to show binding. But a linear 

combination of and , viz., 

cov ionic ' ' 

a o a JD a a d b 

(l-vS) (>^(1)0, (2)+ 0^(2) <!>,(1) +(y-S){<) (IX) ( 2 )+<>, (1)0, ( 2 )> 

_ a 2 a 2 a d d 

{(l-ws) (2+2S^) + (li-s) (2+2S^)+ 8(p-s) ys(l-ys)} 

(3.32) 

should show binding effects for higher values of P / since this 

restores the property of exhibiting higher momentum densities 

perpendicular to the bond-direction than along the bond- direction 

For P = S/ the wavef xonction (3.3 2) is the same as the HL wave- 

function (3.6) and then the rromentum density is the same as given 

by (3.10). Forp>S, (3.32) behaves similar to (3.16) and the 

momentum density is similar to (3.18) . In fact for 1 

(3.32) and (3,16) are exactly Identical. Thus the wavefunction 

16 

(3.32) is expected to show binding for p ^S. Since only non- 
negative values of P are permitted/ an additional condition 
(1-PS) ^ 0 can also be imposed. This leads to the result P l(l/s) . 

It is a well-established result that the overlap integral 
between the two wavefunctions and 'f> 2 (r) is the same as 
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that between their Fourier transforms X,(p) and x^Cp). The 
wavef unctions in the ^-space were depicted by Slater^ as vectors 
and an interpretation given for the overlap integral as the cosine 
of the angle between the corresponding vectors. It is equally 
possible to represent the p-sapce wavefunctions as vectors and 
the relationships in terms of angles between them follow those 
given by Slater^ for the r-space wavefunctions. 

SjTK 

III . 5 SUMMARY AND CONCLUSION 

Momentum densities obtained from the HL and corresponding 
VB wavef unction including ionic terms are compared for diatomic 
molecules. In both the cases, one-electron momentum density 
shows a maximum along the direction perpendicular to the bond. 

The phase relationship for the HL wavefvinction given by the 
oscillator-picture of Coulson^ is not exhibited by the wavefunc- 
tion which includes the ionic terms. However, no physical signi- 
ficance should be attached to this phase-relationship. Since the 
ionic terms in the VB theory play the same role as Cl in MO theory, 
the effect of addinq ionic terms on the CP of H2 was examined. 

The effect of polarization on the CP of H2 was also studied. The 
wavefunction which includes ionic terms is found to give better 
agreement with the Cl- theoretical and experimental profiles, 
especially for the lower q”values. The effect of including 
polarization is not so dramatic, but seems to improve the 
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agreement with the Cl-theoretical profile for the intermediate 
q-values. CP calculated from an MO wavefunction constructed from 
floating AO's is also presented for comparison with the VB pro- 
files. It agrees well with the results for the VB wavefxmction 
which incorporates ionic terms. These calculations point towards 
the importance of including the ionic terms in the calculation of 
molecular CP' s from VB wavef unctions. In Section III. 4, the 
covalent and ionic wavefunctions constructed from OAO' s were 
examined in the momentuaj space. Both the wavefunctions lead to 
an electron density identical to that given by the repulsive state 
wavefunction and, for a given value of the momentum, show minimum 
momentum density perpendicular to the bond direction. Hence these 
wavefxinctions are not expected to show binding individually. This 
is in agreement with the earlier discussion of Slater^ that the 
HL wavefunction constructed from OAO' s does not show minimum dia- 
gonal energy. However, a linear combination (3.32) of the covalent 
and ionic wavefunctions constructed from OAO' s, which leads to 
binding, restores the maxima of the momentiam density perpendicular 
to the bond direction. Thus the criterion of maximtam in the 
piomentum density in a direction perpendicular to a bond is seoi 
to be a sound one for assessing the bonding nature of a two-centre 
VB wavefunction. 
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CHAPTER IV 

ELECTRON MOMENTUM DISTRIBUTIONS AND 
COMPTON PROFILES WITH PSGO MODEL* 

IV, 1 INTRODUCTION 

There has been a spurt in the experimental as well as 
theoretical studies of electron moment\m distributions (EMD) and 
Compton profiles (CP) of molecules during the last five years 
(see Epstein and Tanner^ for a review). In the present chapter 
we shall confine our attention to the theoretical calculations 
of EMD' s and CP' s of some small and medium-sized molecular 
systems. 

The first calculation of the shape of the Compton line 

2 

for a molecule was reported by Hicks using a valence-bond wave- 
function^ of ionic-polarization type for H2 molecule. Duncanson 

*The material presented in this chapter has been/to be pubiished; 
(i) S.R. Cadre, R. Ramaswamy and P.T. Narasimhan, Pramana, 8/ 

99 ( 1977 ); 

(ii) S.R. Cadre and P.T. Narasimhan, Chem. Phys . Letters (to be 
published). 
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and Coulson calculated CP's of methane, ethane, ethylene and 
acetylene using Heitler-London wavefiinctions for the correspond- 
ing bonds. However, the wavef unctions employed by them were of 
limited accuracy. Epstein and Lipscomb^ were the first to obtain 
a general algorithm for evaluating EMD' s from HF-SCF wavefunc- 
tions of polyatomic molecules. This algorithm was used to 
evaluate CP's of boron hydrides^ and hydrocarbons.^ In the 
latter case a localized molecular orbital (LMO) approach was 
used. Some more studies using LMO approach have recently been 

7 

reported for a host of molecules including hydrocarbons. The CP 
of water has been calculated using a variety of wavefunctions and 

O 

compared with recent experimental results. The effect of hydro- 
gen bonding on the theoretical CP of water has also been examined 
9 

in detail. CP's of molecules belonging to 18, 24 and 42-electron 

isoelectronic series have been calculated from ab-initio HF wave- 

functions and the results compared with the LMO CP's.^*^ Thus 

with growing interest on theoretical side and the use of LMO' s, 

13 

the suggestion that CP' s of large molecules may be predicted 
from an analysis of smaller nolecules seems to be coming true. 

There has been a revival of interest on the experimental 
side as well. X-ray CP measurements have been carried out on 
molecule. The CP of lithium hydride has been measured with 
X-rays as well, as yrays. CP' s of a host of hydrocarbons 

have been measured with X-rays. The CP of methane has been 

19 

measured recently with high-energy electron impact spectroscopy. 
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where is the orbital exponent and (A^, B^, C^) are the 
Cartesian co-ordinates of the orbital centre, is used to repre- 
sent an electron-pair, the orbital exponent is related to 
the orbital radius by 



i 


The energy is minimized by optimization of nuclear co-ordinates, 
orbital exponents and some or all of the orbital centres. The 
use of only one Gaussian per orbital is a gross oversimplifica- 
tion and limits the accuracy of such calculations severely. The 
energies obtained with this single-Gaussian version of the FSGO 
model are typically 85% of the corresponding Hartree-Fodc (HF) 
energies. However, such a restricted, sub-minimal basis makes 
this model extremely simple. Due to floating of the orbital 
centres, the FSGO wavefunctions satisfy the vi rial theorem very 
closely. The wavefunctions also correspond to the chemist' s 
picture in terms of bond pairs, inner shells and lone pairs. 

Single FSGO wavefunctions do not represent the inner shell 
properly. Since the inner shells govern the long-range behaviour 
of the CP’s, the CP's computed from single Gaussian version of 
the FSGO model are not expected to give good agreement with the 
experimental as well as the HF-SCF theoretical ones. Use of a 
double Gaussian for each electron pair Instead of a single one 
improves the energies markedly, leading typically to 96% of the 
corresponding HF-SCF energies. Models incorporating a larger 



104 


number of Gau-a-sians have also been studied. A model which 

28 

employs two Gaussians for the inner shell has also been studied. 
Hence it was felt that it would be worthwhile to calculate CP' s 
from single as well as double Gaussian versions of the FSGO model 
and compare the results with those obtained f 2 X)m more accurate 
HF-SCF wavefunctions and experiment. 

5 

Following Epstein and Lipscomb we find that the momentum 
2 

density | x (p) I is given by 


I Pij, f.(p) ^1, m, ® 


-ip. r 






.mj-p-p l^^m^-p' p- 


( 4 . 3 ) 


where Pis the density matrix and fj(p) is the Fourier transform 

til. 

of the radial part of the j orbital. As shown in the above 
g -ip.r., 

work/ e can be expanded in terms of spherical Bessel 


functions 


~ip, 


5:jk 


4 TT £ (-i) 


' jl(pr 








Since only s~type Gaussians are involved in the present model/ 
the radial momentum density can be written, for the single FSGO 
case, as 


i(p) * ®i"0p '^®p 


( 4 . 5 ) 


r pj^ fj(p) fj^(p) V 
j rk 


( 4 . 6 ) 



105 


apart from some numerical' constants. Here i (pr., ) is the 

-^o ^ jk 

Bessel function given by 


= sin(prj^)/(prj^) 


(4.7) 


Pjj, is equal to where is the (j,k)th element of the 

inverse overlap matrix. The Fourier transforms f.(p) and fv(p) 

1 K 

are also Gaussians, 


1 3/4 „ 

fj(p) = (____) exp(-pV4aj) (4.8) 

The spherically averaged CP, J(q) , within the impulse approxi- 
29 

mat ion, can be written as 

00 

J(q) = I ^ dp (4.9) 

lq| 

where q is the projection of the initial momentum of the electron 
on lAe scattering vector. 

Directional CP, JCp^)/ is given by 

JCp^) = / P(p3,,Py.Pz>dPxdPy (4.10) 

which, on substituting for p (p) , yields the following expression 
for J(p„) 5 
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= I 


2T. 

1 j a., a 


3/4 


1 J 


exp(-a. .p ' 
ij^z 


exp( ■ 


4a^^. 7 .2 cos (C^jP^) 


(4.11) 


-11 ( a.+ a.) 

where A^j -|A^-Ajl etc and ^ . Thus, due to the use 

of only s~type Gaussians, J(q) as well as JCp^) can be obtained 
quite easily from single-Gaussian version of the FSGO model. 


CP' s can be obtained similarly from wavefunctions employ- 
ing concentric double Gaussian, viz.. 


= ^il expf-a^Jr-^l 2} + C.j expt-a.2|r-R^|2 J 


(4.12) 


by using Eq. (4.6) and (4.9). 


All the computations for the results reported in this 
chapter were carried out on an IBM 7044/1401 system at 
Kanpur. The integration in Eq. (4.9) was done numerically using 
Simpson' s rule. The total computer time for all the calculations 
presented here was approximately 25 minutes on our IBM 7044/1401 
system. In the next section we present the results of our cal- 
culations of CP's of some molecules. Comparison with other HP- 
SCP data as well as experimental results, wherever avail^le, is 
also made. 
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IV. 3 RESULTS AND DISCUSSION 
IV . 3 . A Hydrogen 

of ^2 calculated from FSGO wavefunction employ~ 
ing just one Gaussian centred at the mid-point of the bond.^^ 

This wavefunction leads to an energy value approximately 85 per 
cent of the HF-SCP energy for H 2 . The FSGC'- theoretical profile 
along with the experimental one^'^ and SCF-theoreti cal profile^ ^ 
obtained from Liu’s wavefunction is shown in Fig. 4.1. The FSGO- 
theoretical profile has J(0) value less than the experimental 
one by about 6 per cent. It also dies out faster than the experi 
mental CP for higher values of q. 

IV, 3. B Lithium hydride 

The radial momentum density, I(p), calculated for single- 
Gaussian FSGO wavefunction is shown in Fig, 4.2. This shows 
an interesting feature in that the l(p) , for the rrolecule is very 
nearly equal to the sum of the individual - orbital contributions 
neglecting the non— diagonal contribution to the momentum density. 
In other words# the curve is nearly identical with the solid line 
in Fig. 4.2 which is the sum of the bond-pair and inner-shell 
contributions. There are_ two maxima appearing in this curve 
corresponding to the maxima of the individual contributions from 
the inner-shell and bond pair. The reason for this additivity 
in the case of LiH is obviously that the off-diagonal terms of 
the inverse overlap matrix, T, are small. The CP calculated 
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from this wavef unction along with the experimental one^^ is 
shown in Fig, 4.3. Our FSGO— theoretical J(0) value is higher 
than the experimental one by about 10% and the width ^ of the 
profile# that is the value of q at which J(q) =^J(o), is lower 
than that for the experimental CP. This is a consequence of 
using a wavefunction of relatively poor energy in which the inner 
shell is not well-represented, 

IV.3 .C Water 

The FSGO theoretical profile calculated from the single- 
32 

Gaussian wavefunction reported in the literature along with 

the directional CP^ J(p ), along tlie C -symmetry axis of the 

molecule is shown in Table 4.1. The CP calculated from a double- 

33 

Gaussian FSGO wavefunction and the experimental CP are also pre- 
sented. The single-FSGO CP dies out faster than the experimental 
one for the higher q-values. In this regard# the present results 
bear a similarity to those obtained earlier for the CP of water 

■ Q 

using minimal basis-set. The corresponding directional CP differs 
from the spherically averaged single-Gaussian CP by about 0,3% at 
the peak of the CP. The CP calculated from the double- Gaussian 
version of the FSGO model agrees better with the experimental 
one. This is# of course# expected because the double-Gaussian 
FSGO' s lead to energy which is typically 96% of the correspond- 
ing HF energy. The inner shell is better represented in the 
double Gaussian version# leading to broader tail of CP, 
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1 ^ 31)1 € ^ ■ 1 




Compton Profile of 

Water 


ISI 

JCp^) along the 
C2-axis.Single- 
Gaussian FSGO 

J(q) spherical 
average. Sin gle- 

Gaussian FSGO 

J(q) spherical 
average . Doiibl e- 
Gaussian FSGO 

4* 

J(q) expt. 

o 

• 

o 

3.387 

3.378 

3.545 

3.73 

0.1 

3.777 

3,368 

3.531 

3.71 

0.2 

3.345 

3.338 

3.488 

3.66 

0.3 

3.29 2 

3.286 

3.416 

3.58 

0.4 

3.216 

3.211 

3.314 

3.43 

0.5 

3.114 

3.111 

3.180 

3.27 

0.6 

2.982 

2.986 

3.016 

3.04 

0.7 

2.835 

2.835 

2.824 

2.77 

0.8 

2.660 

2.661 

2.609 

2.50 

0.9 

2.464 

2.467 

2.379 

2.21 

1 ,0 

2.253 

2.256 

2.140 

1.94 

1.2 

1.081 

1.810 

1.670 

1.47 

1.4 

1,371 

1.375 

1.254 

1.01 

1.6 

0.991 

0.994 

0.992 

0.83 • 

1.8 

0.691 

0.692 

0,678 

0.64 

2.0 

0.475 

0.476 

0.508 

0.52 


^ ^ 8 

+ Experimental Y-ray CP quoted by Tanner and Epstein. The experi- 
mental CP data is given here only upto two decimals. 
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IV. 3 -D Methane 

The single-Gaussian profile computed from a wavefunction 
reported in the literature, the various directional profiles 
along with the double Gaussian CP and the experimental are 

presented in Table 4.2. As noted in the case of water, the 
double- Gaussian CP is in better agreement with the experimental 
one. The directional profiles differ from the spherically avera- 
ged CP by about 0.5% at the peak of the CP. 

IV . 3 . E Ethane 

The single-Gaussian wavefunctions for staggered and eclip- 
sed ethane available in the literature^^ were employed to obtain 
EMD' s and CP's. In Pig. 4.4, the calculated I(p) ys p for 
staggered ethane is shown along with the sum of the individual- 
orbital contributions. These two curves do not agree with each 
other. The reason for this discrepancy appears to be -that v ' the 
non-diagonal entries in the inverse overlap matrix T, are nume- 
rically large in this case and hence, the contribution to l(p) 
from the non-diagonal terms (jT^k) in expression (4.6) is quite 
large. 

We present the various CP's for staggered and eclipsed 
ethane in Tables 4.3 and 4.4 respectively. Here the directional 
CP of staggered ethane, J(P 2 )/ along the direction perpendicular 
to and symmetrical w.r.t. the C-C bond differs at the peak of the 
CP by about 1.3% from the corresponding spherically averaged GP. 
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Table 4.4 

Compton Profile of Eclipsed Ethane 


q/P^ 

'T(P 2 ^ perpendi- 
cular to C-C bond. 
Single- Gaussian 

FSGO 

J(q) spherical 
average. Single- 
Gaussian FSGO 

0.0 

8 . 09 9 

8,202 

0.1 

8.055 

8.145 

0.2 

7.922 

7.975 

0.3 

7.690 

7.695 

0.5 

7.350 

7.308 

0.6 

6.897 

6-824 

0.7 

6.338 

6.253 

0.8 

5.690 

5.615 

0.9 

4.984 

4.931 

1.0 

3.547 

3.546 

1.2 

2.309 

2.330 

1 .4 

1.422 

1.448 

1.6 

0.89 2 

0.916 

1.8 

0.618 

0.624 

2 .0 


0.493/ 
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The spherically averaged CP’ s of these two rotational 
isoiners of ethane can be seen to be quite close to each other 
in both the single and double Gaussian frameworks. ThuS/ within 
the i'SGO— model/ it is not possible to distinguish between the 
CP' s of these two rotamers. 

IV . 3 . P Ethylene and acetylene 

The CP's calculated from single^^ and double Gaussian 
33 

wavef unctions are presented in Table 4.5. As in the previous 
caseS/ the agreement between the theory and experiment is improved 
considerably by the use of double Gaussian FSGO wave functions. 

IV , 3 . G Cyclopropane and cyclobutane 

The single-FSGO CP's calculated from wavef unctions repor- 
34 

ted in the literature alongwith the experimental CP s are pre- 
sented in Table 4.6. The agreement between these spherically 
averaged single-Gaussian theoretical CP' s and those obtained from 
the’ data of Eisenberger' and Marra-^ is not good for higher q-values. 

IV. 3. H Ammonia and hydrogen fluoride 

The CP's calculated from double-Gauss ian FSGO wavefunc- 
tions are presented in Table 4.7 alongwith the experimental or 
HF-SCF theoretical profiles. The agreement of FSGO-theoretical 
values is fairly good/ especially in the higher q-region. 

All the CP's presented in Tables 4.1 to 4.7 are normalized 
to the corresponding number of electrons for O ^q ^15 a.u. or 

o <p^ <15. ’ 



Compton Profiles of Ethylene and Acetylene 
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Table 4.6 

Compton Profiles of Cyclopropane and Cyclobutane 



CYCLOPROPANE 

CYCLOBUTANE 

II 


q 

J(q) spherical 
average. Single- 
Gaussian PSGO 

+ 

E-M 

il J(q) spherical 
!I eve rage. Sin gle- 

i! Gaussian FSGO 

If 

+ 

E-M 

0 .0 

10.480 

10.414 

13.891 

13.836 

0.1 

10.408 

10.318 

13.796 

13 .757 

0.2 

10.193 

10.075 

13.515 

13.435 

0.3 

9.844 

9.662 

13.057 

12.882 

0.4 

9.370 

9.115 

12.432 

12.152 

0.5 

8.782 

8.412 

11.659 

11.270 

0 .6 

8.095 

7.693 

10.738 

10.257 

0.7 

7.324 

6.889 

9.752 

9.185 

0.8 

6.494 

6.059 

8.670 

8.079 

0.9 

5.633 

5.247 

7.548 

6.996 

1.0 

4.776 

4.475 

6.427 

5.966 

1.2 

3.218 

3.155 

4.364 

4.207 

1.4 

2.042 

2.198 

2.774 

2.931 

1.6 

1.304 

1.599 

1.761 

2.13 2 

1 .8 

0.913 

1.250 

1.223 

1.668 

2.0 

0.731 

1.038 

0.975 

1.384 


t Constructed from the corresponding C-C, C-H and C^^, contrxbu- 
tions to experimental CP data of Sisenberger and Marra. 
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Table 4.7 

Compton Profiles of NH^ and HF 


q 


If 

If 

11 

HF 

J(q) spherical 
average .Doubl e- 
Gaussian FSGO 

J(q) I'l 

+ " 

Expt. ;; 

n 

J (q) spherical 
average .Double- 
Gauss i an FSGO 

0.0 

4.215 

4.38 

2.960 

0.1 

4.192 

4.34 

2.952 

0.2 

4.121 

4.27 

2.928 

0,3 

4.000 

4.11 

2.887 

0.4 

3.828 

3.88 

2,827 

0.5 

3.604 

3.58 

2.753 

0.6 

3.334 

3.25 

2.659 

0.7 

3.027 

2.88 

2.548 

0.8 

2.697 

2.50 

2.421 

0.9 

2.359 

2.16 

2,279 

1.0 

2,028 

1.83 

1,126 

1.2 

1.438 

1.30 

1.800 

1.4 

0,988 

0,94 

1,473 

1 .6 

0.680 

0.68 


1 .8 

0.485 

0.50 

0.915 

2.0 

' 0.365 

0.40 

0.709 

+ A. 

Lahmem-Bennani/ private 

commiin i c at ion 
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IV. 4 SUMMARY AND CONCLUSIQWR 

The single- Gaussian FSGO CP' s agree feirly well with the 
experimental profiles, except for the higher q-values. A model 
employing a double Gaussian per electron-pair^^ represents the 
inner shell better and improves considerably upon the single— 
Gaussian EMD' s and CP's. The simplicity of the PSGO model is 
an attractive feature. This is an ab- initio model employing, 
in its simplest version, a sub-minimal basis of s-type Gaussians 
with no empirical or semi-empirical parameters. The orbitals 
used here are "local"^"^ and can be interpreted in terms of classi- 
cal chemical pictures such as inner shell, bond pair and lone-pair. 

25 

Christoff ersen has extended the single-Gaussian FSGO model to 
"synthesize" larger molecules from the corresponding molecular 
fragments. Hence the calculation of EMD' s and CP's of large 
molecules of chemical interest such as various amino- acids, acetyl- 
choline etc. can be undertaken. 

However, as noted earlier, the single Gaussian version has 

severe drawbacks in that it does not represent the inner shell 

properly and leads to poor energies. In this model, the lone- 

pair and pi -bond orbital-centres have to be constrained rather 

3 2 

arbitrarily as they tend to collapse into the nucleus. The 
double-Gaussian version does not have this drawback and gives 
better agreement with the experimental CP' s of molecules having 
pi-bonds or lone-pairs. Tt also represents the inner shell 
better and leads to "tails" of CP v^ich are in better agreement 

* 
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with HF-SCF theory and experiment. 

Since the present model is within the orbital approxima- 
tion, electron- electron correlation effects are not represented 

properly. However, these effects are important in the EMD' s and 

35 

CP's of atoms and molecules. In'spite of these drawbacks, the 
FSGO model offers a very simple way of theoretically calcula- 
ting CP* s for fairly large molecules. 
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CHAPTER V 

isiayil- EMPIRICAL WAVEPUNCnONS 
AND COMPTON profiles 

v-.ii introduction 

Ahlenius and Lindner^ have recently studied the valence 
Compton profiles (VCP) of some small and medium-sized molecules 
calculated from various semi- empirical MO methods. In this 
study# they have examined iterative extended Huckel (lEH), 
CNDO/2 and INDO wavef unctions in the momentum space and the 
corresponding VCP' s. They have calculated various -< p*» values 
for small and medium sized systems and discussed briefly the 
question of transferability of bond profiles in the INDO- 
localized picture. They employed the iterative localization 
procedure of Edmiston and Ruedeiberg. Ihis localization pro- 
cedure produces localized molecular orbitals (IMO) which show 
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marked resemblance to the chemist' s picture in terms of bond- - 
pairs, lone-pairs etc. Ahlenius and Lindner applied this pro- 
cedure to INDO wavefunctions and found that the total VCP* s 
before and after localization agree through four decimals. More 
recently, Figeys, Geerlings and Alsenoy^ have studied the prob- 
lem of rotational invariance of the rnomentum density and radial 
momentum density obtained from NDO (Neglect of Differential 
Overlap) approximate wavefionctions. They found out that only 
following two approximations are physically significant after 
the above two requirements are imposed; 

(i) NDO wavefunctions are used and P(d) and l(p) are 
approximated to INDO or GNDO level. 

(ii) Overlap integrals are explicitly takai into account 

* . 4 ' '"'k 

while solving Roothaan SCF equation together mth unappTO^ 

mated expressions. 

Due to explicit treatment of only valmce electrons, using 
minimal basis-set and employing various approximations to mole- 
cular integrals, semi -empirical calculations are computationally 

■m 

easier 'and much less expensive than the ab-initio ones. Yet the 
one-particle density matrix obtained by semi-empirical procedures 
bears close similarity to Hartree-Fock-Roothaan density matrix 
th many cases. Due to these reasons, they are very suitable for 
handling relatively large, and chemically interesting ^molecules 
and have gained popularity in such calculations. 
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In this chapter, we present our VCP results ojbtained for 
some small hydrocarbons from CNIX)/2 wavefunctions and consider 
specifically the question of bond— transferability of CP's for 
some chemically interesting systems such as allene and butadiyne. 


V.,2 CNDO/2 WAVEFUNCTIONS 

In the CNDO (Complete Neglect of Differential Overlap) 

5 

method only valence electrons are considered explicitly, all 
the inner shells being treated as part of an unpolarizable core. 
The MO is then expressed as a linear combination of atomic 

orbitals (ICAO) 

^. = E t. Oi (5.1) 

^ i ^ ^ 

This leads to Roothaan's equations 

Z Cv( = 0 (5.2) 


or, in matrix notation 


PC = ESC 


(5.3) 


111 (5.2) Pyv is the matrix element of the Fock operator P, 
given by : ; ; 


Pyv 


(5.4) 



H being one-electron operator which includes kinetic energy 

operator and the nuclear-electron potential energy operator. 

/y-l -2 V^(r)}(!.^ dr (5.5) 

being the potential for the atom A due to its nucleus and 
inner shells. G contains the electron -electron repulsion terms 
and 




(5.6) 


where ( v'^fXa ) is the two electanon repulsion integral 


{ vv|Xo)= / •6„(1) (>^,(2) (2) dridr., 

V ^ '1 2 ^ ^ ^ 


(5.7) 


and P is the density matrix where 


occ 

P^« = 2 F c. c. 
A0 IX icr 


(5.8) 


An iterative procedure is then followed to obtain self-consis 
tent MO^ s. 

In the CNDO method, the number of electron-electron repul'’ 
Sion integrals is drastically reduced, by Introducing il 
differential overlap (ZDO) approximat 


zero 
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for arbitrary values of dr and 4 and c^tred on different atoms. 

■ ■ 

Some of the important approximations involved in the QTDO method are: 

(i) The overlap integrals are set to zero unless v . 

(ii) { vv|xo) = o iinless u^vandg = X. in which case 

( pv|xa ) = (yv|xx ) =Y,,x , V 

(5.10) 

(iii) Y^x assumed to depend only on the atoms A and 

B on which w and ^ are centred, and not the actual 
type of the orbitals. Thus we have Y , = ^ for all 

yV Ad 

P and X . 

(iv) Integrals (v 1 Vg |v ) where and belong to the same 
atom , A, are put equal to zero unless v = v in which 
case the integral is taken to be identical for all 
valence AO' s on A 

(v\ Vgly ) = ( 5 . 11 ) 

Vg is the potential for the atom 3 due to its nucleus 
and inner shell electrons. 

(v) The off-diagonal core matrix elements between AO's and 
(}» (which are centred on A and B respectively) are appro- 
ximated by 




o 

AB 



(5.12) 


being the overlap integral aid is a parameter 
depending only on the nature.,-^f A and B» 
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The j^ock-operator matrix elemaits for the CNDO/2 version 
are given by 

^ I ^^bb“ \b (5-13a) 

^ ®AB I ''^AB <5. 13b) 

where ly and A y are the ionization potential and the electron 
affinity respectively of the orbital y , is the total charge 
density on atom A and is the ' core charge' of atom A. Using 
(5,13) calculation is carried out till self consistency is 
achieved in the Pock-operator matrix elements. 

1 

Pople and Beveridge have published a standard CHDO/2 and 
INDO computer program. All the calculations reported in this 
chapter were carried out using this program to generate the 
CNDO/2 wavef unctions. 

V.3 COMPTON PROFILES FROM CNDO WAVEFUNCTIONS 

For the calculation of CP' s the r-space wavefunction is 
Fourier transformed to the p-space 

■tent- ■ . 

i ’J'i(r) e"^£*^ar 


<5.14) 
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The ncmentum density |x(pf is given, 8 apart from constants. 


as 


lx (p)l= 2 Z f,.(p) .Yt . , ’ , , 

vA i j 1 3 (0„» .Y. (6 , (5> ) 

^ Y P p l^.rn^ p p 


where P is the density matrix with 


•ip.r. . 


(5.15) 


occ 


V=1 


(5.16) 


In the 2D0 approximation ^ in momentum- space only diagonal matrix 

elements in Eq. (5.15) survive and the radial mcroen turn density is 
given by ' 

Kp) = ffPiilfiip)! 


As pointed out earlier# Figeys et, al. have proved the rotation- 
al invariance of I(p) for CNDO wavef xinctions . Ihis also leads 

to rotationally invariant CP. The CP# J(q)# within the impulse 

9 . . ' 

approximation is given by 


J( q) = I / I(p) /P 

N 


dp 


- (5.18) 


which reduces to 


J(q) * I Pii 


(5.19) 
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on substitution from (5.1?) . Here J^(q) is the CP due to i^^ 
orbxtal. Hence# the CP can be calculated, from wavefunctions 
employing ZDO approximation, with the diagonal elements of the 
density matrix P and the individual orbital CP' 

For the STO basis used in the CNDO/2 method, the momentum 

space wavef unctions for the individual orbitals are given analy“ 
10 

tically. We present in Table 5.1 various expressions needed 
in the CP calculations using STO, 

It may be noted here that the valence Compton profiles 

(VCP) have been extracted^ ^ from the experimsitally measured 

total CP’s by subtracting out the corresponding atomic inner 

12 

shell contributions. In the literature two different sets of 

values have so far been used for this inner-shell contribution 

to theoretically calculated CP’ s. Epstein has sxibtracted out 

the atomic C^g CP^^ whereas Whangbo et al.^^^^"^ have calculated 

the core contribution by using canonical molecular orbitals. 

These two sets differ considerably, e.g., they differ by 12% 

13 

at the peak. As pointed out by Smith et al, this discrepancy 
between CMO -and LMO core profiles is not due to quality of the 
basis— set but to localization itself. localization leads to more 
contracted core and hence to broader profile. Therefore a compa- 
risdn among GNDO/2 and li^O theoretical as well as experimental 
VCP' s which employ different levels of approximation for the 
core contribution, has to be done cautiously. 



Compton Profiles from SIX) 
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'^•4 RESULTS MD DISCUSSION 

In Table 5.2 we present VGP‘ s of methane and ethane cal- 
culated from CNDO/2 wavefunctions alongwith the profiles cal- 
culated from the data of Eisenberger and Marra^^ and experi- 
mental CP data of Klapthor and Lee.^^ Ihe corresponding UviO 
14 

profiles are also given there. In Table 5 ,3 CNDO/2 VCP' s of 

ethylene and acetylene alongwith those calculated using the 

data of Eisenberger and Marra^^ are presented. Ihe correspond— 

14 

ing LMO VCP for ethylene is also given there. As noted 

earlier by Ahlenius and Linc3ner^ the a greemoit between QJDO 

and LMO VCP’s is fairly good eitcept for higher values of q 

where GNDO/2 profiles are broader than the corresponding LMO 

ones* In the case of methane the agreefr^nt between the profile 

11 

obtained from the data of Eisenberger and Marra# the experi- 
15 

mental CP and the CNDO-theoretical profile is fairly good. 
However, as pointed out in Section V.3/ different sets of values 
have been used in the literature®^ for the C^g contribution to 
CP' s in LMO approach. On the other hand, CNDO/2 method treats 
only valence electrons explicitly. Thus an aribiguity is intro- 
duced in the comparison between LMO and CNDO/2 theory as well as 
experiment. 

Ahlenius and Lindner have pointed out that VCP's are not 
sensitive to localization and the VCP's before and after locali- 
zation agree through four decimals. Hence in the present work 
we have considered it unnecessary to proceed with localization. 
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Table 5.3 

Valence Compton Profiles of Ethylene and Acetylene 



ETHYLENE 


tl 

11 

ACETYLiKE 

q 

CNDO/2 


E-M* 

ft 

ft 

CNDO/2 

0,0 

6,962 

6.996 

6,744 


5.563 

0,1 

6.868 

6.920 

6.664 


5.497 

0,2 

6,594 

6.716 

6.456 


5.304 

0.3 

6.171 

6.368 

6.100. 


5.002 

0.4 

5.609 

5.896 

5.648 


4.615 

0.5 

5.043 

5.312 

5.104 


4.17 2 

0.6 

4.424 

4.672 

4.512 


3.700 

0.7 

3.815 

4.004 

3.904 


3.226 

0.8 

3.241 

3.352 

3.304 


2.768 

0.9 

2.719 

2.752 

2.740 


2,343 

1.0 

2.255 

2.216 

2.052 


1.960 

1.2 

1.513 

1.376 

1.17 2 


1.331 

1.4 

0,992 

0.820 

0.860 


0.880 

1.6 

0.641 

0.484 

0,708 


0.572 

1.8 

0.413 

0.292 

0.412 


0,370 

2,0 

0,267 

0.180 

0,296 


0.239 


'M i l, I' . ' I " H. , — — ■ ^ ' .■ ; , ;■ J; 4' . ' ; ' 

"^Constructed from IMO data of Smith and Whangbo. 

* Constructed from bond- contributions to experimoital CP s 
following Eisenberger and Marra.^l 
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From the VCP-S of CH^, C^H^, and 

bona additivity of GP‘ s we may write 


and consideration of 



^ ^^C-H 




= ^c-c-" 

6J^ 

C-H . 



= ^c=c-" 

^'^C-H 



= ''CzC ^ 

^‘^C-H . 

(5.2®) 


We present CNDO/2 VCP' s for these bonds derived from ( 5 . 20 ) in 
Table 5 •4. By employing these bond VCP* S/' VCP' s for some other 


hydrocarbons may be computed, e.g.. 


C 4 Hg ( bu t adi en e ) 

^*^ 0=0 '^C-c ^'^C-H 


'^C 4 H 2 (butadiyne) 

2Jc=c C H 


'^C 3 H 4 ( allene) 

2Jc=;c '^Jc-H* 

(5.21) 


We have also computed the VCP* s directly from CNDO/2 wave- 
functions for these systems. The geometries of butadiene (tra:is) 
and alien e were taken from ref. 16. The bond length data in the 
case of butadiyne was taken from ref. 17. We present VCP* s cal- 
culated for these systems directly as well as from additivity 
procedure (5 , 21 ) in Table 5.5. We nobe that the differoice 
between the directly calculated CP using CNDO/2 and CP obtaiped 
by addition of bond-profiles is approximately 0*07%, 0,03l& and 
0,18% respectively at the peak of the profile in the case of 
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Table 5.4 

CNDO/2 Valence Compton Profile Contributions, 
from Various Bonds* 


q 

C-H 

c-c 

G=C 

c=c 

0.0 

1.223 

1.026 

2.070 

3.117 

0.1 

1 . 204 

1.018 

2.052 

3.089 

0.2 

1.149 

0.993 

1.997 

3.006 

0.3 

1.066 

0.948 

1.909 

2.871 

0.4 

0.962 

0.891 

1.790 

2.690 

0,5 

0.849 

0.821 

1.647 

2.474 

0 *6 

0.734 

0.743 

1.487 

2.232 

0.7 

0.624 

0,660 

1.318 

1.978 

0 .8 

0.5 23 

0.577 

1.150 

1.723 

0.9 

0.433 

0.496 

0.987 

1.478 

1,0 

0.355 

0.421 

0.836 

1.250 

1.2 

0.234 

0.292 

0.579 

0.864 

1.4 

0.151 

0.197 

0.388 

0.578 

1.6 

0.097 

0.129 

0.254 

0.379 

1 .8 

0,062 

0.084 

0.165 

0.246 

2.0 

0.040 

0.055 

0.107 

0.157 

*Vide Eq. 
electrons 

(5,20) of text, 
respectively.-^:; 

CP's nonnalized to 2, 2, 

4 and 6 


Valence Cbmpton Profiles of Butadiene Ctrans), Butadiyne and Allene 
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butadiene, butadiyn e and allene respectively , We also note 
that the directly calculated CP* s predict soinev^at higher^^^ ^ ^ 

< p > values than those obtained from bond profiles. , - This is 
in accordance with the observation that these molecules contain 
at least one C~C bond shorter than the corresponding average 
C—C bonds, giving the molecule extra stability and thus predict- 
ing lower energy than that obtained from the ' synthetic* profiles. . 
On the whole one can see the remarkably good agreement between 
the CNDO/2 results and those obtained using bond- additivity. However, 
whether these profiles in turn will agree with experiment after 

taking proper care of the core contributions . remains to be seen.. 

18 ' ' 

Epstein has recently pointed out that CP me asuronents using 
high-energy electron scattering techniques, pushed to maxi- 

mum accuracy, may be capable of giving J(o) with an error less 
than 0.1%, Very accurate CP measurements can also be \andertaken 
by scattering of synchrotron radiation. This method is capable 
of giving large scattered intensity with a much better signal to 
background ratio than the other methods, it is hoped that the 
validity of bond additivity for CP's as well as semi-empirical 
wavefunctions for calculation of CP* s may be substantiated in a 
large number of molecules in the near" future with the availabi- 
lity of more refined experimental techniques. 

V;.5 SUMMARY 

VCP* s of methane, ethane, ethyl^e and acetylene calculated 

from CNDO/2 wavefunctions have be^ ’Compared with IMO-theoretical 



and experimental CP's. The agreement among them is fairly good. 
However, as noted in Sections V.3 and ¥.4, the VCP‘ s obtained 
by LMO approach are not nnique due to use of differait sets of 
values of "core contributions'’. These CNDO/2 VCP« s were employed 
to obtain C-H, C-C, C=C and CSC bond VCP' s. 
butadiene, allene and butadiyne were "synthesized" from the 
VCP contributions. The results were then compared with directly 
calculated CNDO/2 VCP' s for these molecules and found to agree to 
within 0.2%. It is pointed out that with the possibility of 
increased accuracy in the experimental roeasuronents of CP' s the 
bond additivity and semi-^enpirical approaches may be evaluated 
more thoroughly and assigned their rightful places. 
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CHAPTER VI 

LOCALIZED MOLECULAR ORBITALS AND 
Y-RAY COMPTON PROFILES OF ACETONE 
AND ALLYL ALCOHOL* 


Via INTRODUCTION 

The localized molecular orbital (IMO) approach (discussed 
earlier in Section 1.9) to Compton profiles (CP) makes it possi- 
ble to break-up the theoretical CP of a molecule in terms of 
CP's of components such as inner shells, bond-pairs and lone- 
pairs. These localized component CP's are transferable and can 


be used to estimate theoretically CP's of several other molecu- 


les by involving additivity of such localized components. Epstein 
and Smith and Whangbo^ have obtained IMG contributions to CP*s 


*The experimental measurements presented 
made on y-ray Compton profile set-i^ at 
Division, Bhabha Atomic Research Centre 
author is grateful to Dr* N.S.,;,^stya 
Physics Section, BARC for roaSclng -evai lab 
mental facilities at BARC. du 

Chaddah for making the autiMa-rybonys^^ft 
at BARC and several hfelpful'disms^-®^^ 
to thank Dr, V.C* Sahni for 


in this chapter were 
the Nuclear Physics 
(BARC), Bombay . Lbe 
thy,; Head, Solid State 
le to him the experi- 
^ tao Mr, Praveai 

their CP set-up 
autdior also wishes 
iltbst -in- this work. 
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from inner shells, bond-pairs and lone-pairs of inolecules con- 
taining first-row atoms from Hartree-Pock quality wavef actions. 

On the experimental side, Epstein et al.^ studied three isomers 
of molecular formula , viz., dioxane, n-butyric acid^^^ 

iso-butyric acid by X-ray scattering technique and compared t^ 
experimental GF s with IMO- theoretical profiles. Ifclt et 
have measured the CP s of cyclopropane and propylene followed by 
a comparison with iMO-theory. Chaddah and Sahni^ have carried 
out a similar study on methyl formate aid formic acid, in this 
chapter we present results of our CP measurar^ts on two Isomers, 
viz,, allyl alcohol and acetone. These are then compared with 
the CP’ s constructed from Epstein’ s^ LMO contributions to CP’ s. 

In the following section, we give the experimental details. The 
present work was carried out using the Y -ray CP set-up at the 
nuclear Physics Division of. the Bhabha Atomic Research Centre, 
Bombay, 

VI. 2 EXPERIMENT ^ ' 

The Y—ray Compton profile set’^up used for the measurenents 
was similar to that shown schematically in Pig, 1.4. The measure- 
ments were carried out at room tonperature with 59,54 keV Y-radi— - 
at ion from a ^"^^Am source of strength lOQ MCi, The scattering 
angle was 165 + 2 . 5 ®. The scattered,,^ray photons were detected 
by a Si (Li) solid-state detector. : -The overall instruitental 
resolution was 650 eV (fuU-wiath at half maxinEm) and o» be 
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assiamed to be quite satisfactory for studying the CP’ s of ace- 
tone and allyl alcohol Which are expected to have a width -v 1.5 keV 
based on IMO approach. The output from the detector was fed to 
a preamplifier, followed by a gated biased amplifier and finally 
energy- analyzed by a multichannel analyzer (MCA) with 400 
channels. The overall stability of this detecting svst©n was 
not very good and readings which showed a drift of more than 
one channel per run were discarded. The MCA employed for accu- 
mulating counts was found to be non-linear in energy and had to 
be calibrated to obtain the energy corresponding to a given 
location of MCA. The calibration curve is givai in Fig. 6,1. 

The energy scale was fixed by using the 26.3 keV and 59.54 kev 

241 

Y-ray peaks from the Am source. The region aroxmd 290230 


channels was found to be very nearly linear. The energy corres- 
ponding to 29 2^^ channel was 48.45 keV ( the CP was centred at'«« 
292^"^ Channel) and the slope in this region was 0.112 keV/ 
channel. Samples were kept in glass cups in an air-tight vessel 
covered with a thin mylar foil. Backgroiind runs were recorded 
with the sample removed. However, during the course of the 
measurements it was found that the backgroiand was quite hi^. 


typically -v 40% of the scattered intensity observed with the 
samples in position in the Compton peak region. The counts for 
acetone and allyl alcohol, wepe awnmula^ and 3600 min. 
respectively, giving -v 8000 and 7500 reafJ^tively m the 

Compton peak region. The backgrotmd t^ts^^i^ this region were 


'''5000 and '''4500 respectively* 


^^ati#tlGal imcertainty 
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in the Compton peak region was of the order of 1.6%. The raw 
experimental data are presented in Table 6.1. Ihe instrumental 
resolution function (IRF) was measured with a l\iCi source 

in order to ensure approximately similar countina rate as in the 
actual experimental me asurenent. The IRF is shown in Fig. 6.2. . 

VI. 3 DATA ANAliYSIS AND RESULTS 

As discussed earlier in Section 1. 5. A, two different 
approaches can be followed for data analysis; 

(i) Deconvolution of the experimental CP with the instru- 
mental resolution function (IRF) and comparing the results with 
theory. ■ ; ; 

(ii) Convolution of the theoretical CP with IRF followed 
by comparison with experiment. 

Paatero et al.^ have studied convolution in CP measure- 
ments in detail. The deconvolution apprpacb poses numerical 

problems due to the algorithm used and yet one has to convolute 

" 6 

the theoretical CP with the residual instrumental function 
before comparison with the experiment. The residual instru- 
mental function may be quite broad and may have oscillatory 
character. To avoid these difficulties, we have chosen the 
latter approach of convoluting the theoretical CP with IRF, 
which was used earlier by Chaddah ^d Sahni. Theoretic^ 

CP' s were obtained following the IKO approach of Epstein.^ %us 



Table 6.1 
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Compton Scattering by Acetone and Allyl- 
Alcohol-Raw Experimental Data 


Channel 

No. 

Energy 

(keVl 


ACETONE 


ALLYL 

ADODHOL 


Total 

counts'^ 

(A) 

Back- 

gromd 

(B) 

(A-B) 

Total 

Counts* 

(c) 

Back- 

ground 

(d) 

(c-d) 

1 

2 

3 

4 

5 

6 

7 

R 

276 

46.660 

2548 

847 

1701 

2167 

1 

762 

S 

1405 

277 

46,772 

2698 

894 

1804 

2459 

805 

1654 

278 

46.884 

2862 

878 

1984 

2650 

790 

1860 

279 

46.996 

3103 

927 

2176 

2907 

834 

1073 

280 

47.108 

3477 

954 

25 23 

3 249 

859 

2390 

281 

47,220 

3858 

971 

2887 

3657 

874 

2783 

282 

47.332 

4374 

1043 

3331 

4034 

939 

3095 

283 

47.444 

4937 

1107 

3830 

4778 

996 

3782 

284 

47.556 

5649 

1223 

4426 

5395 

1101 

4194 

285 

47.668 

6532 

1444 

5088 

6047 

1300 

4747 

286 

41,180 

7254 

1541 

5713 

6585 

1387 

5198 

287 

47.892 

8530 

2153 

6377 

7854 

1938 

5926 

288 

48.004 

9741 

2714 

7027 

8801 

2443 

6358 

289 

48.116 

10892 

3317 

7575 

9820 

2985 

6835 

290 

48.228 

12091 

4291 

7800 

11499 

4375 

7124 

291 

48.340 

12899 

4917 

7982 

11742 

4425 

7317 

29 2 

48452 

13115 

5068 

8041 

11980 

4561 

7419 

293 

48.564 

12896 

5164 

7832 

12015 

4738 

7277 

294 

48.676 

12474 

4963 

7511 

11459 

4467 

6992 
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Table 6.1 (contd. ) 

^ __ _ 


295 

48.788 • 

12016 

4674 

296 

48.900 

10809 

3970 

297 

49.012 

9118 

2892 

298 

49.124 

7695 

2194 

299 

49 . 236 

6833 

1954 

300 

49.348 

5898 

1760 

301 

49.460 

5027 

1595 

302 

49.572 

4158 

1180 

303 

49.684 

3547 

962 

304 

49,796 

299 2 

636 

305 

49.908 

2551 

450 

306 

50.020 

2202 

321 

307 

50.132 

1804 

258 

308 

50.244 

1657 

205 


6 7 8 


7342 

10854 

4207 

6647 ■ 

6839 

9708 

3573 

6135 

6236 

8296 

2603 

5693 

5501 

6749 

1975 

4774 

4879 

5809 

1757 

415 2 

4138 

5254 

1584 

3670 

3468 

4356 

1382 

2974 

2978 

3543 

106 2 

2481 

2585 

3142 

866 

2276 

2354 

2756 

572 

2184 

2101 

2316 

405 

1911 

1880 

1974 

289 

1685 

1546 

1705 

232 

1473 

145 2 

1455 

185 

1270 


+ Data accumulated in 4000 minutes 
* Data accumulated in 3600 minutes 


Counts 
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for acetone we have 


12Jc-h * 


while for allyl alcohol 


Acetone 


= 4Jfn + 2J^, + 4J, 

5s °ls 


0 


+ 


Ip 


'^Allyl 

alcohol 


4Jr 


2Jr 


4J, 


"Is 


1 s 


O. + 
IP 


lOJ 


C-H 


2Jc-0 


2J 


0-H 


( 6 . 2 ) 


The LMO-CP* s of acetone and allyl alcohol calculated from 
equations (6.1) and (6.2) enploying Epstein's^ data are given 
in Table 6.2. 


Data analysis was carried out with the formulae given in 
references 7 and 8. The theoretical CP was converted to energy 
scale by the formula 


(l-cos9)/mc^ 

p - me 5 5 ry:,- {6,3) 

(w^^+ a >2 -201^03 2 cose) 

where me = 137 in atomic units, thus obtaining a fxonction J(a» 2 ) . 
Subsequently corrections for ^sorption in the sample, energy 
dependence of relativistic Compton cross section and energy 
dependence of detector efficiency were applied* I(tij 2 ^ was 
obtained after applying all these corrections. 

1 (^ 2 ^ = A( 0)2) * G( . 0(0)2) • J(« 2 ) ^6.4) 
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T^le 6.2 


LMO- 

Theoretica 1* 

Compton Profiles of Acetone and Allyl 

Alcohol 

q 

Acetone 

Allyl 

alcohol 

— n 

:: g 

If 

ft 

Acetone 

Allyl 

alcohol 

0.0 

13.504 

13.440 

1.7 

2.068 

2.064 

0.1 

13.392 

13.332 

1.8 

1.812 

1.804 

0.2 

13.036 

13.013 

1.9 

1.617 

1,611 

0.3 

12.482 

12.474 

2.0 

1.465 

1.459 

0.4 

11.752 

11.750 

2.1 

1.325 

1.319 

0.5 

10.888 

10.902 

2.2 

1 . 208 

1.202 

0.6 

9.929 

9.952 

2.3 

1.13 2 

1.128 

0.7 

8.901 

8.957 

2.4 

1.027 

1.023 

0.8 

7.871 

7.911 

2.5 

0.962 

0.958 

0.9 

6.888 

6.926 

2.6 

0.904 

0,900 

1.0 

5.952 

5.984 

2.7 

0.845 

0.843 

1.1 

5.112 

5.138 

2.8 

0.795 

0.791 

1.2 

4.376 

4.396 

2.9 

0.753 

0,751 

1.3 

3.732 

3,742 

3.0 

0.717 

0.715 

1.4 

3.188 

3.198 




1.5 

2.743 

2.743 




1 ,6 

2.370 

2.366 





* Constructed from li^O CP data of Epstein,^ 
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Here, A (“ 2 ^ correction due to absorption in a flat 

sample 


T -d(u(ti) ) cosec ctt y(w„) cosec g ) 

A(o.j) = ^ 1 2 


tv (»,) + v(»,). 


(6.5) 


and y(“ 2 ^ linear absorption coefficients, d is the 

sample thickness ( 0.4 cm in the present case) and a and & are 

the angles between the sample surface and the incident and 

scattered beams respectively. The linear absorption coeffi- 

9 

cient data were taken from the tables given by Victoreen. 

The energy dependence of the relativistic Compton cross 
section is given by^ 


C<»2) 


= { 


4m„^ IT' 

o 


: .2^ , “i(VpymoC) ^ 


X 


c 

Z O 


2^. 




( 6 . 6 ) 


The data on energy dependence of the detector efficiency was 
kindly provided to us by Chaddah. The detector efficiency at 
various energies was calculated by interpolation of this data. 

l(oi,^) thus obtained by application of the above correc- 
tions via Hq. (6, 4 ) was convoluted after detector efficiency 
corrections with the IRF givori in Fig. 6.2. The IRF presented 



162 


in Fig. 6.2 was measured for the Compton line at 59.54 keV, but 
the same line shape was used to carry out convolution of the 
theoretical CP centred at'^48.4 keV. The convolution in the 
present case is reduced to a matrix multiplication; 



(6.7) 


where the subscripts refer to channel number, f is the theore- 
tical spectrum, g is the IRP and h is the observed spectnom. The 
theoretical CP* s corrected for the energy dependent terms and 
convoluted with IRF were thus obtained. The multiple scatter- 
ing corrections to Compton scattering data we3re discussed 
earlier in Section 1.4. However, we did not apply these multi- 
ple scattering corrections. Accordingly, the presort CP data 
may be uncertain to 3-5% at the peak. Figs. 6.3 and 6,4 show 
the experimental CP' s of acetone and allyl alcohol respectively 
alongwith the theoretical profiles with all corrections stated 
above followed by convolution with IRF. The agreement can be 
seen to be good in both the cases. However, with the accuracy 
avail^le in the present measurements, comparison between the 
experimental CP's of acetone and allyl alcohol is not meaningful. 

VI.4 SUMMARY AND CONCLUSIONS 

Experimental CP' s of acetone and allyl alcohol were 

' 241 

measured by scattering of 59,54 keV Y-rays fjoom a Am source 

and the results compared with the IMO-theoretical profiles after 







Ur.2S 47.78 48.34 48.90 49.46 

. ; ' ■ E(keV) 

' ^ and theoretical profiles for acetone 
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the latter were corrected for energy dependent terms and con- 
voluted with the instrumental resolution fxmction. Multiple 
scattering corrections were not applied. The agreement between 
theory and experiment is good. However, due to several limita- 
tions of the experimental set-up and owing to accumulation of 
relatively smaller n-umber of counts (implying larger statistical 
uncertainties) no meaningful comparison can be carried out 
between the present experimental CP's themselves. Very few 
attempts have so far been made to measure CP* s of isomeric 
compounds and correlate these with some other molecular property. 
However, with recent revival of interest in this area and use of 
very strong annular Y-ray sources^ it may be hoped that many 
more of such studies would be forthcoming. 


t 
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CHAPTER VII 


calculation op atomic and molecular energies 

PROM experimental COMPTON PROPILES* 

VII. 1 INTRODUCTION 

In this chapter we shall examine the theoretical relation- 
ship between Compton profile (CP) data and one of the most 
important quantities of interest to chemists, viz., energy. 
Relationships between the Compton profile, J(q) and < p^> where 
n = O, 1,.... are available in the literature. The test of the 
relationship between J(q) and <p > in particular has been made 
here by fitting the experimental CP data to a linear combination 
of Gaussians thereby facilitating the evaluation of total energy. 


*The material presented in this chapter has been published in 
the form of a Research Note: S.R. Gadre and P.T. Narasimhan, 
Mol. Phys., 31, 1613 (1976). 
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VII. 1. A Relationship Between Energy and Compton Profile 

1 2 3 

Coulson, Benesch and Smith and Epstein have independently 

derived the relationship between the electronic momentum expec- 
tation values <p^> and J(q)/ the Compton profile iCP), of an 
atomic or molecular system. Coulson's^ derivation is as follows; 

Let p be the nramentum of the electron. The expectation value of 

1 2 

kinetic energy in the N-electron system is then <T>=: -^N <p > . 

Applying virial theorem^ viz., ^©tal ~ ~ ^ ^Total 

total energy of the system. Thus 


"Total 


= -^N<p^> = / p^.P (jg) p^ sine^ dBp d^ip dp 


(7.1) 


Using the definition 


Kp) = / P(p) .p^ sine.. d9^d(|) 


P P P 


(7.2) 


expression (7,1) is reduced to 


00 

=TotaX=‘i /pli(p)-ap 


(7.3) 


Since the CP, J(q), within the impulse approximation,^ is given 

by 


J(q) = \ / 

Ul ^ 

the radial momentum density I(p) is 


(7.4) 


Kp) --2p^^ 


(7.5) 


q =p 
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Hence (7.3) is reduced to 


rotal 


dp 

^ q=p 


(7.6) 


Carrying out the integration in (7.6) by parts 


^Total 


p“.J(p) 1 - 3 p‘^.j(p)dp 


(7-7) 


o O 


If the CP dies out faster than q as q ^ the first term 

in (7.7) vanishes and the expression for is reduced to 


=Total = -3 Jq"j(q)dq 


(7.8) 


Coulson also pointed out that if J(q) is found empirically as 
a sum of the individual component J(q)'s (e.g./ according to the 
localized molecular orbital (UNiO) approach^), also 

additive. Thus the additivity of CP* s implies the additivity 
of bond energies. 


VII. l.B Relationship Between < p^ > and Compton Profile 

A more general expression (7.9) for <p^ > was derived 

3 ■ 

independently by Epstein using an approach similar to that 
of Coulson’ s^: 


<p^> = 2(n+ 1) / q^. j(q) dq , n >0 


(7.9) 
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2 

Benesch and Smith also arrived at this equation (Eq. { 1 , 9 )) 
by an identical approach. They also pointed out that this 
equation is valid for 0 ^n _< 4 by considering the asymptotic 
behaviour of the natural orbitals (NO' s) . For s-type NO' s 


CO 

"v? / ^ (pr).r^dr 

S 'TT is -’O 


( 7 . 10 ) 


where X (r) and X (p) are the NO* s in r and p space respectively 
and j^(pr) is the Bessel function 


jj^Cpr) = sin(pr)/pr 


(7.11) 


For large values of p, (7.10) yields 

where 0(p ^) denotes terms of the order of p Since the orbi- 
tals of higher angular momentum give leading contributions 0(p ^) , 
0(p ^) etc. which die out faster# the asymptotic behaviour of 
l(p) (and hence of J(q)) is governed by the s-type NO's. Using 
the cusp condition for s-type NO* s# viz. 


(7.13) 



172 


Eq. ( 7 . 12 ) is simplified to 

^s^P^ "" 2 Xg(0).p“^T 0(p“^) (7,14) 

yielding for large q values. 


J<q) = 

P^(0) .q ® + 0(q ®) 

(7.15) 

where P (o) = E X y 
0 s ^ 


(7.16) 


being the occupation numbers of s-type NO's. Since the contri- 
bution to J(q) from s-type MO* s is the slowest-varying term of 
the contributions from various NO* s (viz., s, p, d. ..types), 

J(q) dies out as q — ^ . Thus the sum rule (7.9) is valid 

only for 0<n< 4. 


VII. l.C Physical Significance of <p^> 

Physical significance of various <p^"> values was expli- 
citly pointed out by Epstein^. < p*^> is simply the nomalization 
constant. Epstein has made a speculation that.it may be possi- 
ble to correlates p >, the average magnitude of the current 

■ ' 2 
density, to the shielding in nuclear magnetic resonance. <p > 

is related to the electronic energy (vide Eq. (7 .8) ) . Relativis- 

. ' . 4 ■ ' 

tic corrections are proportional . to < p >. itoother observable 
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directly available from the CP is <p ^ J(0). Thus Eq. 

(7.8) and the more general Eq. (7.9) have potential applica- 
tions in chemistry. Energies are conventionally derived from 
spectroscopic and thermodynamic data. However, Eq. (7.S) gives 
a new and novel method of obtaining directly from the 

experimental CP. 


VII. l.D Errors in the Evaluation of <p^> from Comipton Profiles 
.3 

Epstein has examined this problem of evaluating various 
<p^> values from the CP's in more detail. He has considered the 
effects of truncation of momentiom- range and experimental errors. 
The former was studied by him by introducing partial expectation 
values 

Ph 

p^(p^) = 2(n+ l) / p^ J(p) dp (7.17) 

o 

where p^ is the highest momentum at which the CP is measured. 
p^(pj^) increases monotonically with Pj^ and 

=, P^^Ph^ " < p^> (7.18 

Epstein analyzed the behaviour of p^(pj^) for Is and 2p type 
STO* s for which analytic expressions for J(q) are available. 
Denoting the difference{ <p^> - p^(s) } / p^(s) by (where 

s = Pj^/a, a being the orbital exponent), he studied the 
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systematic truncation error A^. From this analysis, he obtain- 
ed the highest momentum values(upto which the CP measurements 
should be carried out) for ensuring an arbitrary accuracy in 
various p^(p^) expectation values. He also noted that a similar 
analysis using a linear combination of atomic orbitals should 
provide reasonable limits of the relative truncation errors and 
the momentum range required for obtaining a given accuracy in 
P^'CPj,). 

3 

Epstein also analyzed the effects of random experimental 
error associated with each point of the CP. He concluded that 
the problem of truncation error can only be resolved by increas- 
ing the range of momentum in CP measurements. With Y -ray scatter- 
ing techniques the present limit is 30 a.u. However, given a 
fixed number of analyzing channels, increasing Pj^ will increase 
the energy spacing between two successive channels and lead to 
greater statistical errors. Ihe statistical errors can be reduced by 
either increasing the counts taken at each point or by decreas- 
ing the spacing between two successive points. Another approach 
would be to record more accurate data at higher momentum than at 
low values of momentum. These points should be borne in mind 
when one makes use of experimental CP data available in the 
literature for evaluating < p^>. 
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VII . 2 FITTING OF EXPERIMENTAL DATA TO A LINEAR COMBINATION 
OF GAUSS IANS 

In the light of the above discussion, it was felt that 
one way to overcome the truncation problen in calculating < 
values from the experimental CF* s would be to fit the available 
experimental CP data to a prescribed analytical form. Another 
approach^ is to append the theoretically calculated high-momen- 
tum atomic "tails" to the experimiental data. 

We- have studied the former approach and fitted the experi- 
mental CP' s to a linear combination of Gaussians: 

J(q) = Z C. exp(-“.q^) (7.19) 

i ^ ^ 

Optimized C^ and can be obtained by a least-squares fitting 
which minimizes the error in fitting given by 


Error = 2:(Jg^p^^(q) - J ( q) ) 


(7.20) 


for all the q-values for which J(q) is available experimentally. 
The advantage of such a Gaussian-fit is that all< p^ values 
(vide Eq. (7.9)) can be evaluated analytically. The energy. 


Ejotal given by 


■'Total 


_ Mu 

4 


C, 

y i_ 

i a 3/2 


(7.21) 
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and 


/ J ( q) dq = 


h 

2 


C. 

1 


a. 

1 


1/2 


(7.22) 


We have examined the CP* s obtained by X-ray, Y-ray as well 
as electron scattering for some atomic and molecular systems. 
With 6-8 Gaussians (i.e., 12-16 parameters), reasonably good 
fits were obtained in all the cases examined. The error in 
fitting expressed by equation (7.20) was typically '\'10 Our 
Gaussian fits obeyed the normalization condition 


/ J(q) dq = N/2 , (7.23) 

o 

N being the nToraber of electrons in the system, to within 0.8%. 

7 

Subroutine STEPIT was employed to obtain these least-squares 

fits. We have also checked our fitting procedure by employing 

8 2 

the theoretical CP for Be 2s orbital. Our< p > value thus 
obtained agrees with the theoretic si value to within 0,2%. 

In the next section we give some relevant details about 
Subroutine STEPIT used for optimization. 

VII. 3 SUBROUTINE STEPIT 

STEPIT is a general optimization subroutine. It can be 
used to find local minima of a real function of several varia~ 
bles. STEPIT uses only function values - no derivatives. It 
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can be employed for optimizing fxinctions which are completely 
non-linear in the parameters- Only restriction is that the 
function to be minimized must be smooth, i.e., it must be a 
continuous function defined on a continuous domain and all the 
first partial derivatives must be continuous. A direct search 
procedure is employed in this routine, v/hich is as follows: The 
direction of the valley is determined by a cyclic relaxation- 
first parameter is varied till the function value is minimized 
with respect to it. The same procedure is then applied successi- 
vely to all the other parameters. When an improvement in the 
function value is obtained, step sizes are increased or decreas- 
ed accordingly. When a minimum is bracketed,parabolic interpola- 
tion is used. After two such cycles of variation of parameters, 
a comparison is made between the direction of resultants. This 
variation of the direct search procedure has an advantage that 
the function value always decreases and the search never diverges 
except in the cases where there is a minimum at «( . To find a 
global minimum of functions of more than one variable, there is 
no known method. A judicious guess of the starting values 

of parameters followed by direct search appears to be the only 
solution to this problem. 

call statement is CAli STEPIT (PUNK) . be 

declared EXTERNAL, is the name of the f unct ion- ccmiputing sub- 
routine, The variables used in the COMMON statement are as 
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In the next section^ we present our results of Gaussian fits to 
experimental CP data. 

VII. 4 RESULTS AND DISCUSSION 

The Gaussian-fit parameters and the corresponding errors 
in fitting are given in Table 7.1. Table 7.2 presents our 
energy results calculated from the experimental CP* S/ the 
corresponding Hartree-Fock (HF) and experimental or the best 
theoretical energies along with the value of the normalization 
integral. As mentioned above, our Gaussian fits obey Eq.(7.23) 
to within 0.8?4. The experimental profiles are generally normal- 
ized by extrapolation to large q-values. However, J(q), in most 
of the small atomic or molecular systems, dies out at q values 
from 5 to 15 a. u. 

The energy values presented in Table 7.2 are given only 
upto two decimals in view of the limited accuracy of the experi- 
mental CP data. We note here that the calculated energies are 
close to but generally greater than the corresponding HF energies 
The deviation from HF energies is upto 6 per cent and can be 
attributed to large experimental errors {of the order of 10-20 
per cent) as well as to the lack of data in the higher q^ regions. 

Klapthor and Lee^^ have recently fitted their methane CP 


to a form 



Table 7.1 
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Gaussian-Fit Parameters (vide Eq. (7.19) of 
text) to Experimental Compton Profiles;^ 


System 



Parameters 



Er 2 X)r in 
fitting^ 

X 10 

He 

C. 

-0.02064 

0.37980 

0.42539 

0.17924 

0.02455 


(X-ray) 

1 

a. 

1 

1.63710 

1.05000 

1.92300 

0.44444 

0.38016 



c. 

1 

0.07101 

0.01200 






a, 

1 

0.19777 

20.00000 




1.22 

He 

c. 

-O:. 23984 

0.31340 

0.85605 

0.04080 

-0.05 285 


( electron) 

1 








a. 

X 

2.09000 

0.43000 

1.72000 

8.00000 

0 . 20000 



c. 

1 

0.11269 

0.03793 






ot. 

1 

0 . 50000 

0 . 1 2000 




0.91 

• Ne 

c. 

0.97174 

0.80807 

0.61127 

0.15129 

0.21545 


(X-ray) 

1 

a. 

1 

0.82665 

0.61200 

0.25389 

0.15761 

0,10526 



c. 

1 

-0.05596 

0.06835 

-0.20383 





a , 

1 

0.56382 

0.04089 

3 .75640 



3.70 

^2 * 

C. 

X 

1.07 230 

0.80091 

0.24686 

0.06224 

0.05435 


(Y-ray) 

ai 

0.61382 

1.58850 

1.08010 

0.38280 

0.12983 



C. 

X 

0.10335 

0.07966 

-0.01815 





ai 

0.14336 

0.77843 

10.00000 



1.78 

Ho 

C. 

0.5 2056 

0.61567 

0.35799 

0.02095 

0.00002 


2 

1 







(X-ray) 

a, 

X 

3 . 297 20 

1.11340 

2.90770 

0.24759 

0.10000 



C. 

X 

-0,00022 







ai 

2.00000 





2.12 

S 

C/ 

1 

0.98892 

G. 25790 

0.89640 

0.91695 

-0.67547 


(X-ray) 

a . 

X 

0.49281 

0.11049 

0.85 299 

2.18040 

3.15920 



C. 

■X 

0.02099 

0.02781 






a . 

i 

2.93450 

3.95550 




3.87 


7* Compton profile data were taken from the following: P. Eisenberger, 

Phys. Rev., h2, 1678 (1970) (He (liquid) : X-ray, H 2 : X-ray) ; P.Eisenberger, 
Phys.Rev.Jfi, 628(1972) (He: Ag X-ray/ O: average of Mo and Ag X-rays) ; 
P, Eisenberger and W. A, Reed# Phys# Rev## 3 237 (1975/ Y -iray) i 

H.F* and R.A. Bonham^ Phys* Rev., A 7^ 1568 (1973) {Hes 

electron). 

± Vide equation (7. 20) of text. ^ ^ 
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Table 1.2 


Energies (in atomic units) from Experimental Com.pton 
Profiles along with the Normalization Integral 


System 

-E( calculated 
from experimental 

Compton profiles; 

(in a . u . ) 

^H-F 

(in a.u.) 

-Et-s the best 
theoretical or 
experimental 
(in a.u.) 

'-k 

J J ( q) dg 
° (total) 

He 

2.69 (X-ray) 

2.87 + 

2.90^ 

1.003 


2.73 (electron) 



1.001 

Ne 

125.16 (X-ray) 

I28.55:j: 

129.94+ 

5.009 

«2 

1.13 (x-ray) 

1.13® 

1.17^ 

0.995 

^2 

105.23 (Y-ray) 

109.00® 

109.56^ 

7.052 

°2 

146.07 (X-ray) 

149.67® 

150.41^ 

7.973 


^ For the sake of \inifo3anity. Is contributions were subtracted out 

wherever necessary from the literature experimental Cbmpton 

12 '' 

profiles. - <P > core electrons from ref . 8 were added up as 

inner shell energies to thos« obtained from the fits. 

* Normalization of the valence profile plus 0.5 per core electron. 

t E. Clementi, IBM J. Res. Dev., 9, 2 (1965). 

a A.C. Hurley in "Advances in Quantum Chemistry", _7, 322, Ed. P.-O. 
Lbwdin, Academic Press (1973) . 

b R. Daudel "Electronic Structure of Molecules", pp. 15, 21, 42, 
Pergamon Press (1966) , 
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(7.24) 


A similar form was also suggested by Epstein^^ as it yields 
correct asymptotic behaviour expressed by Eq. (7.15). This 
form has a distinct advantage over the presently used linear 
combination of Gaussians and should be studied further. 


VII. 5 SUMMARY ME) CONCLUSIQIIS 

Experimental CP* s have been fitted to a linear combination 
of Gaussians. This procedure has the advantage that all < p^ > 
values can be evaluated analytically from the Gaussian-fit para- 
meters. Ihe energies for He, Ne, N 2 / H 2 and 0^ have been evalu- 
ated from X-ray, Y-ray as well as electron scattering CP data. 
These energies are generally larger than the corresponding 
Hartree-Fock energies and this deviation is upto 6 per cent. 

This can be attributed to truncation of CP data as well as other 
experimental errors. It, therefore, seems to us that the energ- 
ies extracted by the use of Eq. (7.8) from the corresponding 

12 

experimental CP's to assess the role of correlation in momentum 
distributions must await further refinement in the quality of 
experimental data (see also Epstein ' ). With the recent revival 

of interest in CP* s and advent of new techniques such as hf-ray 
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scattering, high energy electron impact spectroscopy'' and use 

14 

of synchrotron radiation it appears likely that highly accu- 
rate experimental CP data would be forthcoming. With highly 
accurate CP data and with a better form of fit such as expre- 
ssion (7.24) one may have a novel method of obtaining molecular 
energies which are otherwise usually extracted from thermodynamic 
and spectroscopic data. 
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CHAPTER VIII 


ENERGY AMD COMPTON PROFILE; 
EMPIRICAL CORRELATION IN 
ISO- ELECTRONIC SERIES 


VIII. 1 INTRODUCTION 

There have been very few attempts to correlate molecular 
Compton profile (CP) parameters, such as width of the CP, ^ 
and the peak value J(0) and various <p^> expectation values 
with other molecular properties- Epstein^ has recently specu- 
lated that it may be possible to correlate < p> , the average 
magnitude of current density, available from the CP data to 
some other molecular property such as shielding in nuclear mag- 
netic resonance. However, no such correlation has as yet been 

2 

established. Nhangbo et al. have attempted to correlate guali 
tatively the "molecular size" to the J(0) value in the tai- 
elect3ron isoelectronic series Ne, HF, H 2 O, NH^ and CH^. They 
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3 

employed the definition/ introduced by Pobb et al. of the size/ 

of a localized molecular orbital (LMO) electran 6 as the 
a a 

second rroment measured at the centre of the charge distribution, 
Ra , 


Q 

a 



Ra 


where 




( 8 . 1 ) 

( 8 . 2 ) 


Whangbo et al. calculated the sizes of the above molecules 
belonging to the ten-electron series. They argued that as the 
number of atoms increases in this series/ the charge distribu- 
tion becomes more expanded and leads to sharpening of the CP. 
However/ they did not attempt to correlate quantitatively the 
J( 0 ) value with the molecular size in this iso-electronic series. 
We may note here that this parameter, viz., molecular size is 

purely theoretical, defined within an LMO approach and not 

2 

experimentally measurable. Whangbo et al. have carried out 

the calculation of molecular sizes with a basis-set of double- 

zeta quality. Hovjever, the effects of correlation are not 

included as these calculations are within the Hartree-Pock 

framework. These correlation effects abe, however, isflected 

directly in experimental measurements. Iheir importance in 

4 

the CP calculations has been emphasized. 

Chaddah and Sahni have recently measured the CP* s of two 
isomeric compounds of molecular formula C2H^02/ viz., acetic 
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acid and methyl formate. Ihieir' experimental CP of acetic acid 
was broader (as measured by than that of methyl formate 

outside the experimental error. They attempted to explain this 
observation heuristically in two ways; 

(a) Diamagnetic molecular susceptibility of methyl formate 

9 

is higher than that of acetic acid, implying a larger < r > value 
for the former compound and hence larger <p >and broader CP for 
the acid. 

(b) The sum of bond energies is larger for acetic acid 

than that for methyl formate, explaining the broader CP for the 

acid. 

However, regarding their argumoats, three points should 
be borne in mind; 

( i) These criteria (a) and (b) are only qualitative. 

(ii) They have not been studied for a more general case 
of iso- electronic series of molecules. 

(iii) Many cases of pairs of isomers are encountered where 

the above criteria (a) and (b) are not mutually consis- 
tent, ie., if (a) predicts a broader profile for an isomer A than 

that of B, the cirterion Cb) suggests just the opposite. For 
example, in the case of diethyl ether and n-butanol, criterion 
(a) above suggests a broader profile for the former molecule 
whereas exactly opposite trend is predicted by criterion (b) . 
Similarly in the case of hexene and cyclohexane a broader 
profile is predicted for the former compound by criterion (a) 
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above/ whereas criterion (b) predicts a broader prof iie for the 
latter compound. It is therefore not possible to decide upon 
which of the criteria (a) or (b) is applicable a priori . 

3 

Robb et al. have pointed out that the size and directional 
property of electronic charge distribution must be correlated 
with the total energy of the molecular syst^. In the li#it of 
the results of Wnangbo et al. discussed above, it may thus be 
possible to correlate energy of a molecule in an iso**electronic 
series to its J(0) value. In the following sections/ we present 
our attempt to correlate CP in iso-electronic series with mole- 
cular energy. 

VIJI.2 COMPTON PROFILE PARM1ETERS JCO) AND ^ AND JENERGY IN 
ISO- ELECTRONIC SERIES 

'■ S ' '■ 

We have recently attempted the correlation between energy 
E and <l(0) using theoretical CPVs calculated for CH^, NH^, H 2 O 
and HF. Floating spherical Gaussian orbital (FSGO) wavef unc- 
tions with two Gaussians to represent an electron pair (double 

.7 

Gaussian FSGO model were employed. In Fig. 8.1 we present the 
graph of J(o) ys IeJ within the dotible-Gaussian FSGO model- The 
curve is a smooth one and it will be noticed that J(o) decreases 
monotonically with increasing }Ej . In Fig. 8.2, our graph of 
Jq g , the width of the CP, ys Je} for these molecules within the 
double-Gaussi an FSGO model is presQited. Ihe plot turns out 
to be a straight line in this case ^d Jq ^ 
cally with |e (. 


increases mono ton i- 
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We may attempt to imderstand the reason for the mono tonic 
nature of the curves in the following way. It is known that 
<P is related to the peak height of the profile viz., j{0), 
through the relation 

< P~^ > = y J(0) (8.3) 

If the wavefunct ions satisfy the virial theorem 

- ^ <p^> = E (8.4) 

then we see that in an iso-electronic series, niimerically higher 

2 

energy values correspond to higher < p > . However, with higher 

2 -1 
<p > values, one may qualitatively expect lower <p >values and 

consequently lower J(o) . This trend of higher {e) values corres- 
ponding to lower J(o) values is clearly seen in Fig, 8,1. In 
the ten-electron iso-electronic series considered here, Jq ^ 
corresponds to q values around 1 a.u. This region of q is signi- 
ficant from the energy point of view. It may be recalled (see 
chapter VIl) that the molecular energy E is related to J(q) by 
the relation^ 

“ 9 

E = -3 /q J(q) dq (8.5) 

From the above relation see that the smaller q-value region 
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is not significant for energy. Further at very large q values^ 
we have J(q) dying out as q and hence this region is also 

not siginificant for energy. Hence# it is the intennediate q- 
region which is relevant for energy and the values of q corres- 
ponding to Jq ^ belong to this region. The monotonic increase 
of |e j with Jq g is to be xmderstood in this context. 

We have also tested the snooth correlations between jE| and 

J(0) as well as IeI and Jq ^ on the ten-electron series of atoms 

■ ' 2 “^ +'■ ' ■— ■ ' ■ 2 '-'' 

and ions, viz.. Mg , Na , Ne, P and O . The free ion CP's 
were computed from Hartree-Fock wavefunctions obtained by 
Paschalis and Weiss. For Ne, CP calculated from Hartree-Pock 
wavefunctions^^ was employed. Figures 8.3 and 8.4 show the 
graphs of Hartree-Fock 1 e| ya JCo) and Hartree-Fock | E | vs Jq ^ 
respectively for this iso-electronic series. The general trends 
in these graphs are very similar to those observed above in the 
case of dovible-Gaussian FSGO-CP* s. 

We have attempted to examine the more general empirical 
correlation of J(q) at a given value of q with jE j for a few iso- 
electronic series employing a wide variety of wavefimctions. We 
present the results of this attempted correlation in the follow- 
ing section. 

VIII. 3 CORRELATION OF J(a) WITH MOLECULAR ENERGY IN ISO- ELECTRONIC 

; : 

In Fig. 8,5 we present the graph of theoretical | E j vs J(q) 
at different values of q for ^ calculated within 




Bm0 of HF’-theoreticai /£/ J(0)(in a.uj for tO -electron serie 

or atoms and ions 









197 


the double Gaussian FSGO model. We again observe that all the 
points fall on a smooth curve for a given value of q. ffowever, 
as q increases/ the monotonic nature of the curves changes. For 
lower values of q, the J ( q) decreases monotonically 
increasing jE| . For intermediate q-values, a family of curves 
is obtained in which the monotonic nature changes gradually, 

(to monotonically increasing behaviour). For higher values of 
q, J(q) exhibits monotonically increasing behaviour when i s| 
increases. 

In going down the series HF -> H^O -> -> CH^ , we find 

that the nuclear charge of the heavier atom decreases and lone 
pairs on the heavier atom are successively replaced by bond 
pairs. Since the dominant contribution to J(q) comes from more 
diffuse valence electrons/ it is expected that J(q) at small 
q-values should increase in the series HF — > H 2 O -**■ -> GH^ . 

This is clearly seen from Pig. 8.5. Since these CP' s are nor- 
malized to ten electrons each/ it may then be expected that (in 
order to satisfy the normalization condition) J(q) should 
decrease in the series HF — > H 2 O — > NH^ CH^ for higher values 

of q. 

The behaviour of CP at higher q-values is governed by the 
core electrons. Hiese electrons contribute most significantly 
to the energy. The core energies for these molecules can be 
measured directly by Electron Spect3x>scopy for Chemical .Analysis 
( SSCA) . Since the long-range behaviour of the CP is governed 
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by the more tightly bound inner shells, it may also be possible 
to correlate the SSGA core binding energies. Eg with J(q) for 
higher values of g. In Pig, 8.6 we present a graph of ESCa 
core binding energies. Eg ys experimental J(q) at q= 5.0 for 
CH 4 , NHg., H 2 O and Ne. Ihese points also fall on a smooth curve 
which increases monotonically with Eg. Coming now to the total 
molecular energies in this series, one would expect the core 
energy contribution to dominate. In that event the molecular 
[eI ys J(q) curve should increase monotonically for higher values 
of q. This is also clearly borne out by graphs presented in Pig. 
8,5. Thus it xTiay be expected from the above qualitative argu- 
msits that in ai isoelectronic series J(q) should increase mono- 
tonieally with 1 e I for higher values of q, \diereas it shoxild 
decrease monotonically for lower q-values. For intermediate 
q-region, a gradual transition between these two limiting behav- 
iours is expected. 

In Pig. 8.7 we present a family of curves of jE j ys J(q) 

for various q-values for the lO electron series of atoms and 

2+ + — ' ' ■ -2- ' • t ■ ' . ■ ' ' ■ 

ions, viz.. Mg , Na , Ne, F and O . Hartree-Fock wavefunctions 

and energies were used to obtain the CP and ] E {data. In Fig. 8.8 
a similar family of curves for the 18-electron series consisting 
of C 2 

employed theoretical CP's and energies calculated by Hirst and 

13 ^ ■' ' ‘14 ' , 

Lietomann with waveftmctions ^Ich use polarization. Another 

set of curves is shown for molecules belonging to the 24-electron 


H. 


[g, CHgNH 2 , CHgOH and CH^F is presented. Here we have 








various 
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series using the CP data and energies of Hirst and Liebmann^^ 
in Fig. 8. 9. It may be pointed out here that the energy values 
of some of the molecules for which CP data are available in 
reference 15 lie veary close to each other. Hence we have selec- 
ted four molecules/ viz., CH 2 =CH-CH 2 , CH 2 -CH=NH, CH 2 -GH =0 and 

'■ 9 '''' ■' . 

H-C-NH 2 , which do not have too close energy values, to facilitate 
the clarity in plotting the graph. Ihe graphs in Fig. 8.7 through 
8.9 have the same general features as those observed in Pig. 8,5 
in that I E I vs J(g) points lie on a smooth curve for any given 
value of q. For lower values of q, J(q) decreases monotonically 
with increasing |e | while for higher values of q, J(q) increases 
monotonically with 1 e[ . 

Encouraged by the above observations with theoretical 
energy and CP data employing a wide variety of wave functions, 
we have attempted to use the smooth nature of ) E j vs J(q) 
curves to predict empirically the CP of irolecules belonging to 
an iso-electronic series. Experimental energy and CP data of a 
few other molecules in the iso-electronic series are used first 
to establish the nature of the smooth curve. Then using the 
experimoital energy for the molecule of interest, we cai inter- 
polate and obtain the J(q) values. Since a family of curves! E j 
vs J(q) is available, the entire CP for the molecule of interest 
can be constructed. 

In the literature, unfortunately, experimental CP data for 
molecules belonging to an iso-electronic series are scanty. In 
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f act we were able to locate in the literature, the experimental 
CP data for only one iso-electronic series. This series is the 
10— electron series CH^, H 2 O, HP and Ke. Even here experi- 

data on He are not available. The CP data on CH , 

■ ■ . . ‘ 4 ' ' 

' ' '17 ■' 18 ' ' 

H 2 O and Ne were taken from the literature. We have employed 

in the present vork the experimental CP data of only three mole- 
cules (CH^, H^O and Ne) to establish the nature of | e j vs J(q) 
curves- Prom these curves we were able to emplfically predict 
CP* s of NH^ and HP. In Pig. 8,10 we present the graph of the 
absolute value of experimental energy,} e | ys J(q) for various 
values of q for this iso-electronic series. For each value of 
q, three points corresponding to CH^, H^Q and Ne were connected 
by a smooth curve and experimental energies for NI^ and HP were 
then employed to empirically prelict their CP* s- In Table 8,1 
we present the results o-f the empirical prediction of CP* s of 
NH^ and HP along with experimental and/or Cl-theoretical data 
for these molecules. We note here that in the case of NH^, the 
results are in excellent agreement with the Cl-theoretical CP of 
Ahlenius and Lindner- They also agree, within the limits of 
experimental error, with the experimental CP of Nf^ measured by 
scattering of 35 keV electrons. It has been pointed out to us 
that more precise CP measuremoats are being carried out on NH 3 - 
It is hoped that the empirically predicted CP curve for NP^ 
would be confirmed with higher accuracy. As mentioned earlier, 
no theoretical CP data are available in the case of HP. Our 
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predicted CP agrees well with the CP constructed from the IMO- 
theoretical data of Whangbo et al.^ 

If the presently observed smooth correlation between the 
CP and energy in an iso-electronic series of molecules turns out 
to be e general one, it has yet another potential application. 

The smoothness of J(q) vs jE[ can be used to establish the 
nature of these curves from data on three or four molecules 
belonging to an iso-electronic series. The CP of any other 
molecule of interest belonging to the same iso-electronic series 
can then be measured and tte family of J(q) vs | e{ curves for the 
iso-electronic series can be used to empirically predict the 
energy of this molecule. As the % error in J(q) is rrdniraxm 
around q= 0 , it would be better to obtain energy eiipirically 
from the J(0) vs |b | curve. Thus, the energy of a molecule may 
be predicted empirically, apart from the conventional spectro- 
scopic and thermochsmical methods. 

VIII. 4 SUMMARY AND CONCLUSIONS 

A remarkable empirical correlation that the variation of 
J(q) with (e) is a smooth one for any givai value of q in an 
iso- electronic series has been established. For lower values 
of q, J(q) decreases monotonically with increasing |Ej, whereas 
for higher q-values J{q) increases monotonically with increasing 
j e| . This correlation has been verified to exist for 
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theoretically computed CP* s and the corresponding theoretical 
energies of a host of wavefunctions as well as experimental J(q) 
and experimental {e| . This result has been used to empirically 
predict the CP’ s of NH^ and HP from their known experimental 
energy data, ibe CP' s thxis obtained agree well with the availa- 
ble Cl-theoretical data on these two molecules. The predicted 
CP results of Ni^ also agree with the experimental data within 
the limits of experimental error. Confirmation of our predic- 
tions with regard to molecules in other iso-electronic series 
must await further experimental work. It may thus be worthwhile 
to measure CP's of a number of molecules belonging to an iso- 
electronic series with the same experimental set-up. If the 
empirical correlation illustrated here is verified to be a gene- 
ral one/ as we suspect/ it will have a useful predictive value. 
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